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Abstract 

The bulk S-Matrix can be given a non-perturbative definition in terms of the flat 
space limit of AdS/CFT. We show that the unitarity of the S-Matrix, ie the optical 
theorem, can be derived by studying the behavior of the OPE and the conformal block 
decomposition in the fiat space limit. When applied to perturbation theory in AdS, this 
gives a holographic derivation of the cutting rules for Feynman diagrams. 

To demonstrate these facts we introduce some new techniques for the analysis of 
conformal field theories. Chief among these is a method for conglomerating local primary 
operators Oi and O2 to extract the contribution of an individual primary Oa/ in their 
OPE. This provides a method for isolating the contribution of specific conformal blocks 
which we use to prove an important relation between certain conformal block coefficients 
and anomalous dimensions. These techniques make essential use of the simplifications 
that occur when GET correlators are expressed in terms of a Mellin amplitude. 



1 Introduction 



Exact theories of quantum gravity should be formulated in terms of gauge invariant observables 
associated to the boundary of spacetime. In fiat spacetime, the only such observable is the S- 
Matrix, so a theory of quantum gravity in fiat space will be a theory that computes scattering 
amplitudes holographically. Since AdS / CFT [H El [3] provides a non-perturbative description 
of AdS theories via a dual CFT, one can obtain the bulk S-Matrix from a fiat space limit of 
AdS. This defines a holographic theory for fiat space using a sequence of CFTs with increasing 
central charge. The introduction of the Mellin amplitude ^ [5] for CFT correlation functions 
has led to progress [HI El El El CD] along these lines, and in particular, we recently argued [S] 
that bulk; locality can be understood by showing how the meromorphy of the Mellin amplitud^ 
leads to an analytic S-Matrix. The purpose of the present work is to demonstrate how the 
unitarity of the S-Matrix can be derived directly from the unitarity of the CFT. Specifically, 
we will derive the usual optical theorem 

-i{T-T^)=T^T (1) 

and cutting rules for the 2-to-2 scattering amplitude of massless scalars at a non-perturbative 
level from the conformal block decomposition and the operator product expansion of the CFT. 

Before we outline the derivation, let us first comment on how the standard, manifestly 
unitary definition of the S-Matrix can be applied in AdS/CFT. The S-Matrix is usually defined 
as the overlap between in and out states 

Sal3 = (ain|/3out) (2) 

where a and /3 are multi-particle states composed of asymptotically well-separated, exactly 
stable particles. From this point of view, unitarity arises as a consequence of the completeness 
of the in and out bases, and all of the structure of scattering is encoded in the fact that these 
bases are different. In the interaction picture we write the S-Matrix as 

= («free |5| /3free) whcrC 5 = T 1 6* /"""c" ^^W'^* } (3) 

and T is the time ordering symbol. The unitarity of the S-Matrix follows automatically from 
the unitarity of the S operator. 

All of these statements have simple analogs when we take the fiat space limit of AdS / CFT. 
The key is to realize that global AdS behaves likes a cavity or finite sized 'box' [H], so to 
obtain the S-Matrix we need only setup the correct experiment and then take the size of the 
box to infinity. As originally discussed in [121 US] cind recently revisited in [HI [151 E], one 
can setup initial states corresponding to incoming particles by acting with CFT operators, 
and then measure the outgoing particles with final state operators after exactly one scattering 
event has occurred. Since we want to scatter finite energy particles in the vanishing curvature 
limit of AdS, we need to study bulk states with energy E so that ER — )■ oo as the AdS length 
scale R ^ oo. An elementary feature of the AdS/CFT correspondence is that AdS global time 



^We reviewed aspects of CFT physics and the MeUin amphtude in our recent companion paper [9], and we 
discussed them in detail in [7], so we urge interested readers to consult these references for a review. 
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Figure 1: This figure shows how the AdS^+i cyhnder in global coordinates corresponds to the 
CFTd in radial quantization. The time translation operator in the bulk of AdS is the dilatation 
operator in the CFT, so energies in AdS correspond to dimensions in the CFT. A scattering 
process in the bulk can be set up by acting with smeared CFT operators at an initial and 
final time that are separated by ttR. In the large N limit, a product of n single-trace CFT 
operators creates an n-particle scattering state in the bulk. 



corresponds to radial quantization 'time' in the CFT, as pictured in figure [TJ This means that 
time translations in the bulk of AdS are generated by the dilatation operator D in the CFT, 
so bulk scattering amplitudes involve CFT states of dimension very large compared to 1, but 
very small compared to the central charge. 

In other words, to compute scattering amplitudes using AdS/CFT we setup an in-state 
by smearing with CFT operators at an initial "dilatation time", we evolve the state with 
the dilatation operator D for a time ttR, and then we measure the result at a final time. 
Now it is easy to imitate the usual interaction picture. When studying CFT operators and 
states with dimension small compared to the central charge N'^, we can separate the dilatation 
operator into D = Dq + j^Dj. Bulk perturbation theory and bulk scattering amplitudes can 
be computed using equation ^ with 

Sap = lim (afroe \Sr\ (3iree) where Sr = T I exp 
where now the states a and (3 are created by products of single trace operators, as discussed in 



Di{t)dt 



(4) 
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p^flSl lQ]. From this point of view, the unitarity of the S-Matrix is a direct consequence of the 
unitarity of the CFT. This description of scattering also immediately explains the S-Matrix 
results of [16], namely that to first order in perturbation theory, the bulk S-Matrix is just the 
matrix of anomalous dimensions (a|D/|/3). 

While this procedure looks familiar from the point of view of the bulk, the setup of equation 
(|4]) does not appear very natural in the CFT, nor is it convenient to use for computations. For- 
tunately, in [9] we proved a conjecture of Penedones jB] that gives an extremely simple formula 
for the S-Matrix written directly in terms of the Mellin amplitude for CFT correlators. This 
formula also leads to a nearly trivial relationship between the conformal block decomposition 
of a CFT correlator and the S-Matrix in the flat space limit. Let us now briefly review the 
conformal block decomposition, which can be viewed as a consequence of unitarity in the CFT. 

In any theory whatsoever, one can insert the operator 1 as a sum over states giving 

M^i) = Y,{0,{xr)02{x2)\a){a\0^{x;)0^{x,)) (5) 

a 

in the case of a 4-pt correlation function. In theories with symmetry one can make further 
progress by organizing the states |a) into irreducible representations of the symmetry group. 
In flat spacetime, this means that one can use Poincare invariance to break the sum into states 
of definite energy, invariant mass, and angular momentum and then integrate over the overall 
momentum of the state. In a CFT, we can organize the states \a) of definite dimension and 
angular momentum into primaries and descendants, where the descendant states can all be 
represented via actions of the translation operator on a primary state. If we organize the 
sum in equation ^ so that all descendants are grouped together with their defining primary, 
we have the conformal block decomposition [171 [HI [El 120] 

A(a:.) = 5^i^A/5i(x,) (6) 

of the CFT correlator, where the Pa/ are fixed numerical coefficients encoding dynamical 
information about the theory. The conformal blocks B^^{xi) are the universal functions that 
represent the contribution of a given primary and its descendants to the 4-pt correlator; these 
functions also depend on the dimensions Aj of the external operators Oj, and were recently 
given in Mellin space in [H El [9] for CFTs of arbitrary spacetime dimension. 

The conformal block decomposition can also be viewed as a consequence of the operator 
product expansion, and this is where its power lies. The OPE says that we can write the 
product of two operators as a sum 

where the universal 3-pt function is fixed by conformal symmetry. If we use the OPE 
twice inside a 4-pt correlator, then we can express that correlator as a sum over CFT 2-pt 
functions with coefficients c^^c^^^. But since local CFT operators are isomorphic to CFT 
states, this can also be viewed as a sum over all the states in the theory, as in the conformal 
block decomposition. In other words, the OPE implies that 

PA,i = c'l,c'l, (8) 
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Figure 2: This figure shows how one can conglomerate k — 2 CFT operators in an k-pt corre- 
lation function to obtain a 3-pt function, and then use these 3-pt functions to determine some 
contributions to the conformal block decomposition of a 4-point correlator. This procedure 
makes it possible to use one order in perturbation theory to say something about the next; 
it is precisely analogous to the way that the optical theorem permits the calculation of the 
imaginary part of the S-Matrix using a phase space integral over the product of lower point 
scattering amplitudes. 



We have derived the well-known fact that the 3-pt correlators in a CFT in principle determine 
all the n-pt correlation functions in the theory. 

The rest of this paper will be concerned with making these ideas computationally useful 
and relating them to the S-Matrix. The key to putting the OPE to work is pictured in figure 
|2j we will refer to the process depicted in the bubble at the top of this figure as conglomerating 
operators Oa and Oh into a double trace operator [OaOh\n,i- This makes it possible to use fc-pt 
correlation functions to determine lower point correlators involving multi-trace operators. In 
particular, we can use information about the correlators at one order in perturbation theory 
to compute terms in the conformal block decomposition at the next order, as pictured in the 
second step of figure [2j 

So how do we implement this conglomeration procedure? Naively, one might proceed 
by defining the double trace primary operator as a linear combination of terms of the very 



4 




Figure 3: This figure indicates how the sum over /c-trace operators with dimension A turns 
into a phase space integral over fc-particle states with center of mass energy A/i? in the fiat 
spacetime limit of AdS/CFT. 



schematic form 8^018^02- Then, by imposing that the special conformal generator annihilates 
the sum, one finds relations for the coefficients. By itself, this is a rather involved combina- 
torially exercise; some partial results were obtained in [211 E], and for completeness we give 
a recursion relation for the coefficients in the case of a general double trace primary operator 
in appendix [C} However, it turns out that determining these coefficients is actually the easy 
part, because to use these coefficients to compute correlators involving a double trace primary 
[Oi02]n/ "we also need to differentiate a CFT correlator involving Oi and O2 a total of 2n + £ 
times. This procedure is very cumbersome, especially at large n and i. 

Fortunately there is a better method that exploits the simple properties of the Mellin 
representation for CFT correlators. Instead of differentiating single trace operators, we can 
integrate them against simple 'wavef unctions' that conglomerate the single-trace operators into 
the desired double-trace state. To form an operator of dimension A and spin £ from two single 
trace operators inside a correlator, we write 

{OaA^)---) = j d''xid^X2fAA^,^1^^2){OiiXi)02ix2)...) (9) 

where '. . . ' indicate any other local operators that may appear in the correlator. The wave- 
function /a/, which we will determine in section [2| has only power law dependence on the 
differences between the coordinates. Because of the simplicity of f^/, when we represent CFT 
correlators as Mellin amplitudes, the integrals in equation ^ can be done immediately using 
the Symanzik star formula, which one can view as the Mellin-space analog of the formula for 
the Fourier transform of e*^'^. We will also see how to use these methods to extract the coeffi- 
cient of an individual conformal block from the Mellin amplitude. In section [2] we will derive 
these techniques and use them to obtain some new results about CFTs, and then in section |3] 
we will also make essential use of this technology in our derivation of unitarity. It seems likely 
that these techniques can be usefully applied far afield from our discussion of the fiat space 
limit of AdS/CFT. 

The process of conglomerating operators at one order in perturbation theory and then 
combining the results to give information about the next order should remind the reader of 
the way that the optical theorem 

-i{T-T^)=T^T (10) 
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computes the imaginary part of the S-Matrix. In section |3] we will show that in fact, the 
conformal block decomposition as computed along the lines of figure [2] reduces to the imaginary 
part of the S-Matrix in the flat space limit of the bulk AdS theory dual to the CFT. 

To derive the optical theorem, we need to show that the sum over fc-trace operators in 
the conformal block decomposition reduces to a phase space integral over fc-particle states, 
as pictured in figure [3] for k = 2. We explain this essentially kinematical fact in section 3^ 



One might also wonder whether all operators that can be exchanged in the conformal block 
decomposition are really fc-trace operators, and what role is played by the operators dual 
to unstable particles. The S-Matrix connects in and out states composed of exactly stable 
particles, and so scattering amplitudes between unstable particles are not well-defined and do 
not appear in the unitarity relation. The qualitative difference between stable and unstable 
particles emerges only in the flat space limit, when the original primary operators dual to 
unstable particles get lost on the sea of multi-trace operators with which they mixes. We saw 
an explicit example of this phenomenon in [9j, where we obtained a Breit-Wigner resonance 
from the fiat space limit of AdS /CFT. 

This also means that the small black holes that can occur as intermediate configurations 
in scattering processes are not literally states in the theory, since they too are unstable. Thus 
there are no 'small black hole operators' being exchanged in the conformal block decomposition, 
and we are not missing any contributions when we formulate unitarity purely in terms of stable 
multi-particle states. 

The conformal block decomposition provides an expression for the exact 4-point correlator, 
but the left hand side of the optical theorem only involves the imaginary part of the S-Matrix. 
Another way of saying this is that in general, the optical theorem does not provide sufficient 
information to fully determine the next order in perturbation theory, because the real part of 
the S-Matrix cannot be uniquely computed. But this means that when we use the OPE as 
pictured in figure [2| we must be missing terms that correspond to the real part of the S-Matrix! 
The missing terms can be most easily understood by looking at the one-loop example in figure 
|4} In the conformal block decomposition of the 4-point correlator computed by this loop 
diagram, when double-trace operators are exchanged there are terms that correspond to 'cuts 
at the edge of the diagram'. These combine the 3-point functions of mean field theory (ie the 
CFT correlators that follow from a free theory in AdS) with interacting 3-point functions, as 
pictured on the left and right sides of figure |4j We will prove that in the fiat space limit, these 
terms in the conformal block decomposition only contribute to the real part of the S-Matrix, 
and so they drop out of the optical theorem. These edge cuts identically represent the part of 
the S-Matrix that is non-trivial to reproduce using dispersion relations. 

Interestingly, the proof that these edge cuts only contribute to the real part of the S-Matrix 
requires an identity first conjectured in [22], which amounts to the statement that the edge 
cut terms are total derivatives. Alternatively, the conjecture says that the OPE coefficients 
Cn,£ = Cn/ + Scn/ foT these edge cuts satisfy 

1 d 

Cn,e5cn,e = --^{cl i-f{n,i)) , (11) 

where 'y{n,i) is the anomalous dimension of the double trace operator [Oi02]n/ of dimension 
Ai -|- A2 + 2n -|- £ + 7(n, £), and Cni and 6cn£ are respectively the infinite N value and finite 



6 




Figure 4: This figure provides a schematic depiction of how a 1-loop Witten diagram in AdS 
decomposes via the conformal block decomposition in the dual CFT. For illustrative purposes, 
the bulk theory has both a a-^id a f x^"^^ interaction. The dashed lines indicate 'cuts'; the 

central cut, highlighted in purple, provides the familiar imaginary contribution to the optical 
theorem in the flat space limit. The conformal block decomposition also includes the 'edge 
cuts' on the left and right, which have no analog in discussions of the cutting rules. These edge 
cuts are very important in order to obtain the full correlator, but in the fiat space limit they 
only contribute to the real part of the S-Matrix, and so they drop out of the optical theorem. 



corrections to Cn/- We will precisely state and prove this statement and a relevant gen- 
eralization at a non-perturbative level in section 2^ by using the conglomeration techniques 
discussed above. 

The outline of the paper is as follows. In section [2] we derive our technique for conglomer- 
ating operators and apply it to some useful examples, obtaining a few new results along the 
way, including the infinite conformal block coefficients in arbitrary spacetime dimensions 
and the generalized derivative relation indicated in equation (11). In section [s] we show how 
the unitarity of the S-Matrix follows from the conformal block decomposition, as we briefly 
outlined above, and we give a full one-loop example. Finally in section |4] we conclude with a 
discussion of the implications and opportunities for further work. 
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2 Conglomerating Operators 



The goal of this section will be to understand how to 'conglomerate operators' in order to 
combine a pair of local primary operators into a third composite operator that appears in the 
OPE of the first two. We will also be able to use these techniques to extract specific terms 
from the conformal block expansion. 

2.1 Basics 

Before we see how the conglomeration process works, let us first review a few basic facts 



about CFT operators. We have written out the full conformal algebra in equation ( 133 ) in the 
appendix, but for our present purposes it will suffice to consider the commutation relations of 
the dilatation operator the momentum generator P^, and the special conformal generator 
K^. These take the form 

[D,PJ = P,„ [D,K;\ = ~K^, [P^,K,] = -2{T]^,D + tM^,) (12) 

The crucial feature to notice is that P'^ and K^^ act as raising and lowering operators with 
respect to the dimension, which is the eigenvalue of D. Unitary CFTs have lower bounds on 
the allowed dimensions of operators, so this means that after some number of applications of 
any state of definite dimension will be annihilated. A state that is annihilated by K^^ is 
called a primary state, and the operator that creates this state is referred to as a primary 
operator. All states of definite dimension and angular momentum can be classified as either 
primaries or descendants of a primary. Since the momentum P'^ = id^, descendants are just 
derivatives of primaries. 

Let us begin with some very concrete examples in mean field theory (ie the dual of a 
free theory in AdS), where all correlators are determined by the 2-pt functions of single-trace 
primaries, which are the operators dual to the fields in AdS. Given two single-trace primary 
scalar operators Oi and O2, one can form double-trace primaries [Oi02]n,e which will have 
dimension Ai + A2 + 2n + £ and spin i. The most trivial example is O1O2, which is primary 
and has n = i = 0. But we can also form the operator 

A2{^^0^)02-A^O^{^,02) (13) 

One can check using the conformal algebra that this operator is primary. A slightly more 
complicated example is the operator given by the linear combination 

- + (m) 

which is also primary, and has n = l,i = 0. In appendix[C]we present a recursion relation that 
completely determines the appropriate coefficients for any double-trace primary with arbitrary 
n and i. At large n, the double trace primary operators approach a one-to-one correspondence 
with the space of 2-particle states in c? -|- 1 dimensions, a fact that will be important later on. 

A very natural question follows: given an n-pt CFT correlator involving Oi{xi) and 02(2^2); 
how do we extract an (n — l)-pt correlator with the double trace primary [Oi02]n,i{xo)7 We 
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could proceed by using derivatives as above, but this quickly becomes extremely cumbersome. 
Furthermore, when we move beyond the mean field theory limit, all operators will pick up 
anomalous dimensions, and when these become large it is difficult to define precisely which 
operator we intend when we write [Oi02]n/{xo). So instead of using derivatives, let us try to 
use an integral over a wavef unction f\ £, and define 

[C»i02]n,£(xo) = j d%d%fA,+A,+2n+eAxo,yi,y2)Oiiyi)02iy2) (15) 

This applies to the mean field theory case, but for general CFTs we can use dimensions A 
other than Ai + A2 + 2n + i. Now we need to determine the wavefunction /a/. Fortunately, 
this can be easily accomplished with the introduction of so-called shadow operators. For any 
primary scalar operator O of dimension A, we wish to find a shadow operator O of dimension 
d — A so that 

{Oix)diy))=6\x-y) (16) 

If we simply define 

(17) 

then one can check that the desired identity is satisfied. Now we can use these shadow operatos 
to compute the wavefunction /a,^. As a warm-up, assume that i = 0. If we compute the 



correlator of both sides of equation (15) with (9i(x 1)02(3^2) then we find that 



/a,.o(xo, x„ X2) oc (^^^)A-a.+a.(^^^)A-a.+a. (18) 

where Xij = Xi — Xj. We did not even need to do any integrals to compute this wavefunction, 
because the 3-pt correlator of scalar primaries is determined uniquely up to an overall constant. 

Before we make use of this result, let us first generalize it to the case where i > 0. For this 
purpose, it will be simpler and more elegant to use the embedding formalism [251 [21], where 
we can use the machinery that was developed and nicely explained in |25]. The basic idea of 
this formalism is extremely simple - since the conformal group in d dimensions is S0{d,2), 
it is most natural to use coordinates that transform in the fundamental representation of 
this group. Thus we will represent each coordinate Xj with a + 2 dimensional vector Pj, 
constrained so that = and identified projectively so that Pi ~ APj for real A > 0. These 
coordinates correspond to the null cone that is the asymptotic limit and boundary of AdS 
when it is regarded as a hyperbola in a + 2 dimensional embedding space. If we use light 
cone coordinates for the Pi and choose the specific normalization P.^ = 1, we find 

{P;',Pr,Pn = il,x^,x^) (19) 

This means that the inner product of the Pi is 

2P, ■ P, = i^+P,r + prp+ - Pl^P^^ = (x, - x,f (20) 

Conformal transformations of the Xi simply act as fundamental S0{d,2) transformations on 
the Pi. We will normalize the 2-pt functions of single-trace scalar primary operators so that 

(0(Pi)0(P2)) = % (21) 
-'12 
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where P12 = 2Pi • P2 and the normahzation 



c r(A) 



with 2h = d, the spacetime dimension of the CFT. As shown in [25], the correlators of operators 
with spin can be described in the embedding formahsm as polynomials in auxiliary d + 2 
dimensional vectors Zj which soak up the tensor indices of the spinning operators. We will 
only make use of the simplest examples from [23], such as the 3-pt function of two scalar 
primaries and a spin i primary 

ZA,...ZA,{0,{P,)02{P2)0iy''^{Ps)) = {c'le)TtA.(Z,P3;Pi,P2) (23) 

wliPrP T^'^ (7P-P P\ - ((Z • Pl)P23 - {Z ■ P2)Pvi)^ 

Wnere J^^^^^l^Z, 1^3, /"l, i"2j = Ai+Az-A+I A2+A-Ai+<; A+Ai-A2+f K'^V 

P 2 p 2 p 2 

-'12 -'23 -'31 

Note that the universal function 7'a/a2 fixed by symmetry, while the 3-pt function coefficient 
provides dynamical information about the theory. The auxiliary coordinates Zj are taken 
to have the property that Zi ■ Pi = 0, and correlators must have a 'gauge invariance' under 
Zi Zi + aPi for any a. One can see immediately that the scalar-scalar-spin-£ correlator 
satisfies this gauge condition. 

We can use these results to determine the general wavefunction /a/. The operator Oa/ is 
defined by 

ZA,...ZAPti'''{P,) = j d''P,d''P2 [ZA,...ZAj^y^'{Po,Pl,P2)] Ol(Pl)02(P2) (25) 

If we again take the correlator of both sides with the product of shadow fields Oi{Pi)02{P2) 
then we find the result 

Zai--- ZaJa/"''^'{Po, -Pi, -P2) = —r'^d-Ai,d-A2i^^ -^o; -Pi, -P2) (26) 



defined in terms of the universal function from equation (23), where A^^^ ^ normalization 
factor that we will determine later. Note that this wavefunction also depends on the spacetime 
dimension d and the dimensions Ai and A2, although we have suppressed this dependence in 
the notation. One could continue on and use the results of [25] to obtain wavefunctions 
involving several operators with spin, but for our purposes equation (26) will be sufficient. 
Now let us see how to use this result to compute interesting 3-pt functions and to extract the 
coefficients in the conformal block decomposition. 



2.2 Using Conglomeration 

The intuition we used above to introduce conglomeration was perturbative. In general, our 
conglomeration procedure can be understood in terms of the operator product expansion. 
Isolating a single term in the OPE, we can write 

0,{P,)02{P2) = ckW^A/(^i' P2)OaAPi) + ■■■ (27) 
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where the '. . .' contain all the other operators in the OPE, including the descendants of Oa/- 
If we compute the correlator of both sides with 0/^/{Po) all of the terms in the ellipsis vanish, 
and we relate the three-point correlator to the OPE coefficient multiplied by the normalization 
of the OPE and the operator Oa/- For example, in the notationally simple case £ = 0, we 



—A + — 

define &a o ~ ^ ^ ^^"^ so we find that 



(Oa,o(Po)0i(Pi)02(P2)) = ci^PE^a^'^" (0A,o(Po)0A,o(Pi)). (28) 
When Oaa is normalized to give a two-point function Pqi^; then we have 



"A,£ ~ '^A/;OPE (29) 



using the definition of the 3-point correlator in equation (23 ). Conglomeration makes it possible 
to extract both 3-point correlators and OPE coefficients. 

Before proceeding to calculate we need to set the normalization. We will use a convention 
such that the general 2-point correlator of Oa/ with itself is |25j 



((Zi ■ Z2){Pi ■ P2) - {Z2 ■ Pi){Zi ■ P2)) 



{OaAZu Pi)OaAZ2, P2)) = ^a:;a.;2;+2. (30) 



-'12 



where we have effectively defined the operator OA,e by conglomeration in equation (25). Now 
we can determine the normalization of the wavefunctions iV^ ^ in equation ( 26 ) by demanding 



{OA/OA,i} = J d'^PifAAPa, Pu P2)fAAPb; ^3, A) (Ol(Pl)02(P2)Cl(P3)C2(P4)) , (31) 

where we have suppressed the dependence of /a/ on the auxiliary variables Zi for notational 
simplicity. One could also set the normalizations in terms of the conformal block decomposi- 
tion. This follows because 

J d''p,d''P2fAAPo; Pi, P2)BA'APi) « - A') V- (32) 

If the result is non-vanishing it can be used to fix the normalization of the wavefunctions; this 
also means that we can use conglomeration to uniquely identify terms in the conformal block 
decomposition. 

Finally, let us consider what happens if there is more than one operator with the dimension, 
angular momentum, and global charges of Oa/- If there are many such operators, they can 
certainly mix with each other, so they can only be differentiated based on their correlation 
functions. By applying conglomeration to different correlators, such as 

(01020102), (01020304), (03040304) (33) 

we can extract all the information we need to separate the operator 0a/ that couples to 0i 
and 02 from the operator 0^^ which has a 3-pt function with 03 and 04. This may be 
relevant for more complicated CFTs. 
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Now that our wavefunctions are normahzed, we can compute 3-point correlators via 



j d^Pid^P^fAAPo; Pi, P2) {0,{P,)02{P2)0^{P^)0,{P,)) = {0^,,iPo)OsiP3)0,{P,)) . (34) 

If there is no operator or operators O^^g in the OPE of Oi and O2 then the result will be zero. 
This procedure is most tractable when the correlators are expressed in terms of the Mellin 
amplitude M{Sij) d El El El E] , so that 

/ioo ^ 
[d5]M{5,,)\[T{5.,)P;;^\ (35) 
i<j 

The reason is that the only dependence on the Pij is in the form of a power-law, and the 
projective integrals over the Pi that we find when we conglomerate can be easily accomplished 
using the Symanzik star formula, which states that 

/n „ n 

d'P n r(/0(-2P. ■ P)-'' = vrM [d6] n m^P^' (36) 
i=i i<j 

where the 5ij integration variables are constrained by 'Ylljkj ^ij ~ ^i- Note that when n = 3 this 
means that the Mellin space integration variables 6ij on the right hand side are completely 
fixed, so there are no integrals to do. One can think of this very useful formula as the analog 
of the Fourier transform of e"^^'^ in momentum space. 

We will first make use of this technology to determine the conformal block decomposition 
of mean field theory in any number of dimensions. This simple result has been obtained for 
d = 2 and c? = 4 in |22], but we are not aware of it appearing anywhere in the literature for 
the case of general d. The relevant correlator is simply 

mP,)02{P2)0,{P,)02{P,)) = 1^ (37) 

-MS -'24 

To extract the conformal block decomposition for £ = and A = Ai + A2 + 2n, we simply 
need to integrate 



1 CaiCa.2 

2d-Ai-A2-A A1+A-A2 A+A2-A1 ^ pAi pZ 

P34 ' P04 ' P03 ' 



rid p rid p t^Ait^A2 /oo\ 

Ct J 3U, J 4 2d-Ai-Ao-A A1+A-A9 A+A0-A1 ^ ^A, nAo ■ V / 



We can apply the Symanzik star formula of equation (36) to the integrals over Pi and P2- In 
both cases the constraints completely determine the integrals over the 6ij, and we find 

2,T{-n)T (h - Ai) Tih- A2) r (-/i + n + Ai + A2) ^ P^2 



r(Ai)r(A2)r(/i + n)r(2/i-n-Ai-A2) V ^^Pot"'"^ot^V ' 

where as usual 2h = d. The coefficient outside the parentheses is c^q^ao- 

In the limit that n approaches a non-negative integer, the above expression is singular. 
Ultimately, we are interested in extracting double-trace operators whose dimensions are exactly 
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given by integer n, and thus one might be concerned about whether we are really able to 
regulate this singularity. For the reader who is interested in such subtleties, we will show 
in Appendix |A] the details of how we choose our regulator. For the more casual reader, 
however, the idea of the following derivation is relatively simple: once we look at the physically 
normalized operator, the singularity in the three-point function cancels against a singularity 
in the normalization factor, so that the physical three-point function is finite. In practice in 

the following, this cancellation of singularities will take the form r(— n)/r(0), which we will 

(-1)" 



take to be '•"Y To fix the normalization factor, we can conglomerate again to compute the 
2-pt function of |C'iC2]n,o- 

In fact, it is worthwhile to pause and note that beginning with any 3-pt function we can 
conglomerate Oi and O2 to obtain a 2-pt function. This relates the coefficient c^q ^^at sets 
the size of the 3-point correlation to a coefficient c^'^ in a 2-point correlator. To compute the 
relation, we multiply the 3-pt function by /a/ and integrating over Pi and P2 to obtain the 
2-pt function 

A,o ^ ^ 7r^^r(o)r {h - ^±^) r(A - h)r [h - ^±^) 

Nl, r(/.)r(^±^)r(^±^)r(rf-A) ' ^ ^ 

We see that there is again a singularity of the form 'T(O)". If the correlator that we are 



computing is of the form in equation (31), then we must have c^'*^ = 1, in which case we find 
the very useful fact 



_ ,12 ^"r(o)r {h - r(A - h)T {h - ^±^) 



Note that this is a non-perturbative result, and is not restricted to mean field theory. 

Now we can incorporate the normalization q and compute the desired conformal block 
coefficient for mean field theory. As we discussed near equation (jsj), the coefficient is simply 



(-12 ^2 ^ r r (Al)n(A2)n (1 + Ai - h)^ (1 + A2 - h)^ 

[c^,o) LA.CA. ^, ^^^^ + A2 - 2/. + 1 + n)n (Ai + A2 - /I + n)„ ' ^^^^ 

where we recall that the Pochhammer symbol (a)^ = r(a + 6)/r(a), and the spacetime dimen- 
sion in the CFT is 2h. In appendix [B| we generalize this method to arbitrary spin conformal 
blocks in the scalar four-point function. The integrals can also be performed in this case, 
with just a bit more book-keeping to track the various terms that appear when we expand the 
degree i polynomial in /a,^. This gives the resulting compact form for the conformal block 
coefficients in mean field theory (ie a free scalar theory in AdS, or a CFT at infinite A^): 

/_12n2 ^ CaiCa2(-1)'(Ai -h+ 1)„(A2 -h + l)„(Ai),+^(A2)^+n 

™^ i\n\{i + hUAi + A2 + n-2h+ l)„,(Ai + A2 + 2n + £ - l)KAi + A2 + n + £ - ' 

(43) 

This result matches that of |22j in the cases they considered, namely that of 2-dimensional 
and 4-dimensional CFTs with Ai = A2 normalized without the factor of CaiCaj- 



^Briefly, one can precisely regulate the r(0) and T{—n) singularities by taking the dimensions A, A' of the 
conglomerated operators [OiO2]a,0 and [C3C4]a',o to differ until physical quantities are calculated. 
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2.3 Further Applications 



While extracting the OPE coefficients of double-trace operators in the infinite theory is a 
useful example of conglomeration, its full power comes from the fact that it is an essentially 
non-perturbative technique, and one can use it to extract the coefficient of an arbitrary operator 
in the OPE. Thus, if we know all the four-point functions of a set of operators, Oi, O2, O3, O4, 
then we can conglomerate O1O2 to make an operator [OiO-2\i\/ of arbitrary dimension A and 
spin £, with A a free parameter. We do not need to know A a priori. Rather, the result of 
conglomerating will give vanishing OPE coefficients except at the values of A for which there 
actually is a corresponding operator in the O1O2 OPE. We will now turn to examples where 
we look at the connected four-point functions of the theory and use conglomerating methods 
to extract information about specific conformal blocks. Mellin space is an essential tool in this 
study, since by construction it organizes the connected correlators into contributions of definite 
powers of the Pj/s. Correlators can then be integrated against the wavefunctions simply by 
repeated application of Symanzik's star formula. 

In this subsection, we will first discuss general results on the application of the wavefunc- 
tions to connected four-point functions, in particular how to extract both OPE coefficients 
and anomalous dimensions. We will then turn to the application of these results to specific 
AdS models. A direct consequence of our methods will be the proof of an important derivative 
relation, discovered empirically in [22], between OPE coefficients and anomalous dimensions 
of double-trace conformal blocks: 

^t>S = ^|^(Ky'7H)> (44) 

where c^^^ are the infinite A^ OPE coefficients of double-trace conformal blocks and 5c^(^ are the 
differences between the exact OPE coefficients and the infinite A^ OPE coefficients. In general, 
this formula is true only to leading order in perturbation theory, but as we will explain in 
detail, for a certain class of contributions it actually holds exactly. 



2.3.1 OPE CoefRcients from Connected Diagrams 

Consider the Mellin amplitude for a four-point function ((9i(920304), and let A^ = Ai = A2 
while Afo = A3 = A4. The four-point function has only two independent Mellin variables, 
which we can choose to be 5 = 2(Aa — ^12), 7 = 2(5i4 + 812 — A^) 

= (Ol(Pl)C»2(P2)03(P3)O4(P4)) ^ Mi234(5,7). (45) 

One can think of 5 in analogy with the mandelstam invariant s, and so when we look in the 
s-channel, the angular momentum information will be carried by the 7 variable. Conglom- 
erating O1O2 to produce [CiC2]a,£ involves integrating the correlation function against our 
wavefunction /a,^. Mellin space is ideally suited for this integration, since its form is already 
a decomposition of the correlator into powers of Pjj's, for which the wavefunction integrations 
just involve a repeated use of Symanzik's star formula. Thus, one obtains a general formula for 



the three-point function of [CiC2]a/ with C3, C4, of the form in equation (23) with coefficient 



4^ = / ^Mi234(5,7)i^A,.(5,7), (46) 
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[Oi02]a,i defined this way still has to be normalized, for which one must calculate its two-point 
function as we discussed in the £ = case near equation (41). In general this is singular, and 



must be regulated by instead calculating the two-point function ([CiC2]a,^[C'iC'2]a'/) with 
A 7^ A'. One then takes A — A' at the end of the calculation, and then only in physically 
normalized three-point function coefficients. In order for the result to be non-zero, we must 
have an oo/oo behavior as A approaches its physical value. We will see shortly that this occurs 
only at values of A for which there is a pole in the Mellin integrand, as we should expect on 
the general grounds discussed in [TJ [9] . 

To avoid unwieldy formulae, we will focus on the special case where the connected four- 
point function in question contains only conformal blocks of spin-0, for instance corresponding 
to s-channel scalar exchange in AdS, and leave the general case to appendix |B.2 Then, the 
Mellin amplitude does not depend on 7, and the four-point amplitude takes the form 



dSd-f (-1) M{d)T{Aa - f )r(Afe - l)T' 



• 5+7 ^ 
. 2 ' 



{2m) 



4 



12 



34 



23 



I 2 



(47) 



We can obtain from this the ([Oi02]a,oC^3C'4) three-point function by conglomerating Oi and 
O2 together: 



a,(p: 



1 



N 



fl2 
A,0 



dPidP2 



P\ 



-A^iPi] 



(4J 



12 



0202 
-'01 -'02 



This may be evaluated by applying Symanzik's integral twice, which introduces two new Mellin 
variables (two from the Pi integration, and none from the P2 integration). However, three of 
these integrations are purely kinematic, in that M does not depend on them, and so can be 
done independently of M{6). Fortunately, performing first the dj integration, all three of them 
take the form of Barnes' Lemmas, and can be computed in closed form. We arrive at 



^3(P.) 



2 ' 



1 



X 



d^ 



r A 



4ivS rWr(2/.-A)r(A)plp#p-^#.A. 



r 



34 
5 



h 



A 

2" 



r A 



M{b). (49) 



So we have obtained a simple formula for the OPE coefficients in this special case where 
and the Mellin amplitude is independent of 7 



Ca,0 



— TT 



2h 



r2(/i 



t)r^(f; 



4ivro'r(/.)r(2/.-A)r(A) 



db 



r 



A 

2" 



r 



h 



A 

2" 



r A 



r A 



M{6). 



(50) 



One can obtain the normalization from equation (41), which can be computed term-by- 
term in perturbation theory if such a series is available. vVe give the general formula for these 
OPE coefficients with arbitrary i and a 7-dependent Mellin amplitude in appendix B.2 
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2.3.2 Conformal Block CoefRcients 



Let us discuss an application of this formalism to the extraction of conformal block coefficients 
and anomalous dimensions of double-trace operators at leading order in perturbation theory 
for a simple AdS theory. Such examples were studied in [221 IISl ES], where different methods 
were used. While [HI |26] found a fairly simple method for extracting anomalous dimensions, 
the calculation of conformal block coefficients remained cumbersome, to say the least. We will 
begin with the decomposition of the four-point function into double-trace conformal blocks 
at leading order in although we will see later that this decomposition also has a non- 

perturbative interpretation: 

A = J]PiH5a„(x,) + lpo(n)7H|^i?A„(x.), (51) 

n 

where Poin) = cl^oC^% and Pi{n) = c1^q6c^q are, respectively, the infinite and correction 
terms to the double-trace conformal block coefficients. The partial derivative with respect to n 
brings down logarithms, since i?A„ in position space contains terms with Xj's to the n-th power. 
When we conglomerate the four-point function with A = A„, we pick up the contribution from 
a specific double-trace operator. We will see later that this form is appropriate for the study 
not just of the leading order in perturbation theory, but furthermore for non-perturbative 
corrections to a large, important class of conformal block contributions that we will refer to 
as "cuts through the edge of a diagram" . Thus, the reader should keep in mind that although 
the specific examples we will compute in this section are perturbative, the general formulae 
we obtain will be applicable for gaining non-perturbative information about the CFT. 

The extraction of OPE coefficients in the presence of anomalous dimensions is a bit subtle. 
Let us therefore begin with a conceptually simpler case, where = '^^"^^^ and A;, = 
are unrelated to each other. Then, O1O2 and do not have any of the same double-trace 

operators in their leading order OPE, so no anomalous dimensions will appear in the four-point 



function at this order. From equations (50) and (40), we find the following expression for Pi: 



G{Aa+n) 



-1 r^(A„ + n) 

^^""^ ~ 4f(0) r(2(A, + n))T{2Aa + 2n - h) ^^^^ 
xvr^ f d5 T{l-Aa-n)csc{7r{l-A^))T{l + Aa + n-h) 
J 2ni ^ ' r (f - A„ + 1) r - A, + 1) sin (tt (f - A,)) 

Because of the singular r(0)^^ prefactor, when we perform the 6 integral by contour integration, 
we can discard any residues that are non-singular when n is an integer. This is a great 
simplification, because the only such contributions are those poles at 

- = Aa + m, m = 0,...,n, (53) 

where there is a pole from both the esc term and the first F function in the integrand. These 
residues are in one-to-one correspondence with the residues of the following equivalent, but 
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simpler, integral, where the F ^(0) has been cancelled: 

' 4 sill (ir (A. -At)) J 2irj ^' F (f - A„ + l) F (| - At + l) ^ ' 

For any specific Mellin amplitude of the form M{6), this formula for the OPE coefficients is 
relatively simple to use, since now for any n it is just a finite sum of non-singular residues. 

2.3.3 Anomalous Dimensions and the Derivative Relation 

We will now generalize the results in the previous section to include cases with anomalous 
dimensions. This will requires addressing the subtlety mentioned above. To see the issue 
explicitly, recall that the perturbative three-point function for two single-trace operators Oi, O2 
of dimension Aq with a double-trace operator of dimensions A + 7(A) takes the form 

. C„o + ()C A+^(A) 

•^3 = 'a. ^ ' (55) 
-'12 



^ (cnfl + K% + 4^^7(A) logM + . . . ^ ^ 



P12 

where u = j^^p—- The problem is that the position-dependence of a three-point function from 
a Mellin amplitude is completely fixed, and cannot contain any logarithms. This is because 
at any finite order in perturbation theory, the anomalous dimensions naively appear to break 
the conformal invariance, and it is only the resummation of all order of the logarithms that 
reproduces a conformally invariant correlation function with the new, shifted dimensions of 
operators. 

However, we are actually in a position to get around this difficulty very easily with the use 
of the results we have just obtained above. While there were no anomalous dimensions for 
Aa and A;, unrelated, in the limit of A;, — )■ A^ we should be able to see logarithms reappear. 
The important physical point to note however is that now there is no difference between the 
double-trace operator with dimension 2Aa + 2n and the one with dimension 2Ab + 2n, so the 
physical three-point function will be the sum of both of these Pi(ri)'s. As one can see from 



the csc(7r(Aa — A;,)) prefactor in eq. (54), each of these Pi(n)'s is singular in this limit. But, 
this singularity exactly cancels in their sum, and the subleading (finite) term in the three- 
point function contains a logarithm! Evaluating this explicitly, we find the sum of the two 
three-point functions as Af, Aa is 



1 



liin (piHm^'"+" + (A,^ A,)) 

9 d \ , f ^^^^^^^T{l-Aa-n)T{l + Aa + n-h) 




dAa dA 



J ^ 'J ^ ^ F(f-A, + l)F(f -A, + l) 



APt^'^dnK ^ " 'J ^ ^ r(f-A„ + l)F(f-A, + l) ' ^ ' 
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Comparing this to the expected form of the three-point function in eq. (56), we see that we 
have derived a simple formula for the anomalous dimensions 7(n) and OPE coefficients pi{n) 
when Aq = Afe: 

^ ^ 2 ^ 'J r(f-A, + i)r(f -A5 + 1) ' 

Piin) = l^PoinMn). (57) 

This proves the relation between OPE coefficients and anomalous dimensions that was 
found empirically in [22]. While we have focused in this section on cases with only spin-0 
conformal blocks, the proof in fact generalizes straightforwardly to any spin. The reason is 
that this result depended only on two properties of our expression for the OPE coefficients: 



first, that only the singular residues in eq. (53) survive the P ^(0) prefactor, and second 



that these residues depend only on A^ through the combination A^ + n, except for the factor 
P (I — Aa + 1) P (I — Afo + 1) that is symmetric in (Aq <h- A;,). This allowed us to exchange a 

d 

by the derivative. 



derivative A^ for one in n, since derivatives acting on this symmetric factor are cancelled 



2.3.4 Example Computations 

Let us now apply this formula to some concrete examples. The simplest possible AdS inter- 
action that affects only scalar conformal blocks is a X(j)^ interaction, which corresponds to a 
Mellin amplitude that is just a constant. Let us now apply this formula to some concrete 
examples. The simplest possible AdS interaction that affects only scalar conformal blocks is 
a interaction, which corresponds to a Mellin amplitude that is just a constant. Then, the 
anomalous dimension is simply 

7(n) oc G{Aa + n)Y2-^ \ P(2A, + n + m - /i) 

m=0 ^ / \ / 

1 {h)n{2Aa + n-2h + l)n{2A + 2n)_h , . 

(A<, + n)t,(2A, + n-/.)„ ' ^ > 



where we have used the expression for the infinite OPE coefficients from eq. (42) and 
dropped an overall n-independent prefactor. This quantity was computed using alternate 
methods in (TB], whose results can be seen to agree with that above. 



Next, let us turn to t<;/>^x^) which has different fields on the left and right and is one 



4 



of the contact interactions in figure 
first order in A; these will be usefu 



We will compute the OPE coefficients 5c2^^_,_2„ q 
ingredients when we study an example of the optical 



theorem in section 3.4 In this case, we apply equation (50) to the trivial Mellin amplitude 



M = A4 = A^P(A2 — /i) n^=i r(A-) • '^^^ normalize the wavefunction by using equation 
(41) with the mean field theory c^*^ +2n0' -^^^ normalization of the wavefunction we find 



../xx_ "^'r(o)r(/^-f)^P(A-/.) Ca,(a,).(i + a,-/.), 

P(/i)P(f)'p(2/i- A) ^/n\ (/i)„ (2A^ -2h + l + n)n (2A^ -h + n)^ 
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where A = 2A^ + 2n corresponds to the dimension of [C'x^xln.o to first order in perturbation 
theory. To apply equation (50) we need only integrate using Barnes' Lemma, giving 



6c 



A,0 



A. 



X 



(-l)"+ir(-n + - A^)r^(A^ + n)T{A^ + + n - h)T{2A^ + n - h) 
2n!CA^r(2Ax + 2n)r(A^ + A^ - h)T{2A^ + 2n - h){A^)n{A^ -h + l)n 

n\{h)n{2A^ + n-2h + l)„(2Ax + n-h)n 



(60) 



We can use this coefficient and the equivalent one with (p ^ ip, to compute one-loop conformal 
block coefficients, as pictured in figure |4j This will be useful for verifying the unitarity relation 



in the fiat space limit when we come to section 3.4 



3 S-Matrix Unitarity from CFT Unitarity 

In this section we will derive the optical theorem 

-i{T-T^) = T^T (61) 

for 2-to-2 scattering of massless scalars by analyzing the conformal block decomposition in the 
fiat spacetime limit of the dual AdS theory. The derivation will occur in several steps. First, 
in section 3.1 we review our recent result from [9], where we showed that in the fiat space limit 
of AdS/CFT, conformal blocks correspond to delta functions in the center of mass energy with 
a definite angular momentum i. This means that 



J2PA,iBi{6,,) 5ocP^,,Cf-^)(cose) (62) 

We will carefully compute the phase in the normalization of the blocks in order to precisely 



identify their imaginary parts. Then in section 3.2 we will compute the left-hand side of 



equation (61) in terms of the conformal block decomposition. Terms in the conformal block 
decomposition that we call 'edge cuts' only contribute to the real part of the S-Matrix[^ these 
edge cuts are pictured in figure [5] and also shown in a perturbative example in figure |4j The 
conformal block coefficients are simply products of OPE coefficients, so that schematically 

Pvs,i = {^"/s, + Sc'^s,i) {c%i + Sc%,) (63) 

In this expression, the edge cuts are simply the terms that involve the mean field theory OPE 
coefficients c^,^. To prove that the edge cuts drop out of the optical theorem we will make 
essential use of the relations we derived in section 12.31 for the double-trace conformal block 
coefficients and their anomalous dimensions. 

The central cuts pictured in figure |5] do contribute to the imaginary part of the S-Matrix, 
so it remains to show that these are equal to the right hand side of the optical theorem. The 



We will refer to the left-hand side of equation (61 1 as the imaginary part of the S-Matrix, although in fact 



it can be complex if the in and out states are distinct. 
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central cuts are exactly the conformal block coefficients given by the product of interacting 
OPE coefficients, so we have that 



-^(^-^0U/= E '^^aA^. (64) 

But the right hand side is already in a form that can be interpreted as the right hand side 
of optical theorem, T^T. It only remains to argue that the OPE coefficients are proportional 
to scattering amplitudes, and that the sum over exchanged operators corresponds to a phase 
space integral over multi-particle states. The OPE coefficients Sc^/ can be computed by 
conglomerating k + 2-point correlators into 3-point correlators involving fc-trace operators, 
and these A;-trace operators provide a basis for scattering states in the flat space limit |15j . We 



will show in section 3.3 that the sum over fc-trace operators becomes a phase space integral 



over fc-particle states in the flat space limit of AdS [15] , as depicted in figure 3 so that 



oo „ k ,^ 



-7 ll(2vr)^2E. 



This will complete the derivation of the optical theorem. Finally, in section 3^ we show how 
this logic applies in a complete one- loop example. 

3.1 Conformal Blocks in the Flat Space Limit of AdS/CFT 

The conformal block decomposition will be crucial to our proof of the unitarity of the holo- 
graphic S-matrix. As we discussed in the introduction, a conformal block corresponds to the 
exchange of a particular state in the CFT, and in perturbation theory conformal blocks sum up 
to give the various different cuts of Feynman diagrams. So, we will ffist review the flat-space 
limit of conformal blocks from [9], paying close attention to the phase in the normalization 
of the result. It was shown there that the flat-space limit of a block is proportional to a 5 
function in the center of mass energy. We explicitly performed the transformation between the 
Mellin amplitude and the S-Matrix P, |9] 

1 f 



T(si,) = lim / rfae^a'^-^^M ( 5,,- 



4a 



1=1 i=l 

where As = | Yli ^« ^r half the sum of the external dimensions. Alternately, one can argue 
based on general principles as follows. Conformal blocks are just the contribution to the four- 
point function (or n-point functions, more generally) from complete irreducible representations 
of the conformal group. Consider a conformal block whose primary has dimension A. The 
primaries are the lowest weight states of the representation, i.e. those annihilated by the 
generators of special conformal transformations. There is a one-to-one mapping of AdS 
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states to CFT states, as well as of generators of the conformal algebra to generators of AdS 
isometries, which in the flat space limit is just the Poincare algebra. Furthermore, special 
conformal generators in the fiat-space limit become the momentum operators [161 US] • Thus, 
all states in the conformal block map onto a complete irreducible representation of the Poincare 
group that has center-of-mass energy y/s = A/R. Consequently, in the fiat-space limit, the 
conformal block; can contribute only at this value of s. Its angular dependence is further 
constrained by symmetry to be the appropriate polynomials in cos 6*, which are Legendre 
polynomials in d = 3 and Gegenbauer polynomials C^'^~^\cos6) more generally. 

To read off the normalization, it is simplest to use the inverse of eq. (66 ), because integrating 
over delta functions is very easy. So we compute 

M(5,) '=""^ff ^e~^fi^--'T = -4/3|) (67) 

and input the form for T{sij) required by symmetry: 

T{sij) = AA,f5(s- A Vi?')Cf"') (cos ^). (68) 



The integration over /3 is then trivial to perform: 



A 



T 



5i2)'^-^^e5^Cf"^^(cos^) (69) 



This result should be compared to the large 6ij limit of the conformal blocks themselves. They 
are fixed up to an overall normalization by conformal invariance, and we will normalize them 
in accordance with our definitions from previous sections, so that 

e"(^-^^) (,..(^+.-2.) _ 1) r(^)r(A - + 1) r (^) r 
""^^^''^ = ^^ r(A.-f)r(A,-i)'^^-^^^^-^' 

where we define Ai = A2 = A^ and A3 = A4 = A;,. Here, r = A — £ is the twist of the 
conformal block and Pe^ri^ij) is a Mack polynomial [H [9]. For £ = 0, it is just Pq.a = 1- We 
will need only the large 6ij limit of the block, with 6ij oc Sij: 

Pi,M,) geA-Si2Yci'~'\cose), (71) 

where the proportionality constant gi^ is real and (70, a = 1- Expanding B^{5ij) at large 5ij 
and A, we obtain the approximation 

Ui=i r(Ai + i-h) 

x(^i- cot (^(2^12 - 2Aa + r))) sin(7r5i2) sin(7r((5i2 - A„ + A;,)). (72) 

Because of the sin and cot factors, this does not strictly speaking have a well-defined large ^12 
limit. However, if we smooth over an 0{1) region of ^12, the last line averages out to 

-i cos(7r(A^ - Ab)) + sin(7r(A^ + Afc - r)) ^^^^ 



21 



Re + Im Real 




Figure 5: This figure depicts how the 'edge cuts', which correspond to the terms in the 
conformal block decomposition involving free propagation, only contribute to the real part 
of the bulk S-Matrix, while other operator exchanges contribute to both the real and the 
imaginary pieces of the S-Matrix in the flat space limit of AdS/CFT. 



As an aside that will be relevant shortly, note that for the double-trace operators r = 2Aa + 
2n-\- 1 and r = 2Ab -\-2n + I, this simplifies further to 



- ^e(-i)'^^-(^-^^) and - ^e^-^y'<^^-^^) ^ 



(74) 



respectively. Thus we have obtained the overall normalization coefficient for the flat space 
limit of the conformal blocks 



A^ 



2N As-ft.-l' 



TT 



3/i-l25+2Aj:^2/i-5^2-/i 



nti r(A, + i-h) 



(75) 



X 



-i cos(7r(A„ - Ab)) + sin(7r(A„ + Af, - r)) 



Note that both terms inside the parentheses appear to be even functions when we switch 
Aa -H- Afc. However, in the special case where the conformal block is a double-trace operator, 
it is important to note that when we make this switch we must take r from 2Aa -t- 2n to 



2A{, + 2n, as was already made manifest in equation (74). Thus we observe the crucial fact that 



the imaginary part of the coefficient for double-trace operators remains even under A^ ^ A;,, 
while the real part is odd. 



3.2 The Imaginary Part of the S-Matrix 

Let us use the tools we have developed to derive the optical theorem. As we saw in the previous 
section, when we take the flat space limit of the conformal block decomposition of the 4-pt 
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correlator we find 



Ti2^34(s.,) = 5^ Pa/ AA,iS{s-AyR')cr'\cose) (76) 



The normalization factor was obtained in equation (76), and Pa/ is the full conformal 
block coefficient. 

In general, all operators Oa/ with the charge of O1O2 can contribute to this sum, and so 
for each dimension A we can write 

Pa, = cl,cZ (77) 
We have expressed the conformal block coefficients in terms of CFT 3-pt functions by using 



the OPE, as discussed in section |2.2[ Now we would like to isolate the double-trace operators 
[Oi02]n/ and [O^O^n/- These operators play a special role because they are present even 
when the bulk theory is free, and so they are responsible for the '1' when we write = 1 + iT. 
In the special case of these operators, we can write 

cS = c5 + (78) 

where c^^^ is the 3-pt function coefficient corresponding to a free bulk theory and 5c^^^ is the 
change in this coefficient due to interactions. We will use a similar notation for c^'*^, although 
note that [Oi02\ni and [02,0^n/ will, in general, have different dimensions and so give distinct 
conributions. Now we can write the conformal block coefficients for these operators as 

p _ 7,12 -34 I (-12 r 34 , r 12 -34 \ , r 12 r 34 

If the operators Oi and O2 are the same as and then the first term corresponds exactly 
to the '1' part of the S-Matrix, and otherwise it is absent. The final term comes purely 
from interactions, and actually combines two different pieces, one involving the exchange of 
[Oi02\n,i and the other involving the exchange of [O^Oi\n,- However, the terms in parentheses 
are precisely the edge cuts pictured in figure |5] They combine free propagation on one side 
with interactions on the other, and are associated with poles in the F functions from the Mellin 
integrand rather than poles in the Mellin amplitude itselfj^ 



This definition of "edge cuts" as any term with, a c factor should be intuitively reasonable, but we can also 
more rigorously connect it to bulk diagrammatics. For any diagram, we can formally label all internal field 
lines by (/>j-^'''s, which are distinct from the external fields (jif^^ This is just a relabeling and does not change 
the Mellin amplitude itself. However, it is now manifest that "cuts through the middle of the diagram" are any 
conformal blocks for operators made of internal fields (ti^P , and "cuts through the edge of the diagram" are any 
conformal blocks for operators made of external fields (f>\ . Since (j)] s never appear as internal lines, such a 
diagram is not sensitive to any lower order corrections to their conformal blocks, and must take the form of a 
leading correction, i.e. 

A D cl2,fc5B2A„+2„+5ci2,c5B2A.+2«, (80) 

plus a possible 5 (0^^)^7125^ i32A„+2n term if [Oi02\n,i — [Oz04\n\i- This is a long-winded way of saying that 
in this labeling, it is manifest that "edge cuts" are exactly equivalent to terms that contain c^^ or c'^^. However, 
c's are exactly identified as the parts of the OPE coefficients that are zero-th order in any bulk couplings, and 
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Let us show that these edge cuts drop out of the imaginary part of the S-Matrix. Assume 
for simphcity that we are deahng with real scalar fields, so all bulk couplings are real; the 
generalization to complex couplings is straightforward. Now the OPE coefficients are real, and 
i(T"'' — T) = 2Xm(T). Furthermore, in any physical process in the flat-space limit, there will 
be a finite resolution much greater than the AdS curvature scale, which we can incorporate 
through a resolution function fe{s, sq) narrowly peaked on s = sq. Taking the fiat-space limit 
of the conformal block decomposition and integrating against this resolution, the edge cuts 
contribute as 



oc 



:^/,((2A, + 2n)V/?2 



-M-irA.^C'f "i)(cos^) + (l,2, A„ ^ 3,4, Afe) (81) 



where indicates the finite resolution and we have used our computation of the normalization 
factors Aa/, defining = (A^ — Af,) for convenience. Only the phase of A^^g is relevant here, 
so we have discarded a real overall coefficient. We are primarily interested in the imaginary 
piece of T for unitarity, but it will be enlightening to keep track of both its real and imaginary 



pieces. We will now use our formula from equation (|54|) for the OPE coefficients. This is specific 
to 



0, but the generalization of the following step to non-zero spins is straightforward: 



T. 



£=0 
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+ 5^/,((2A, + 2n)Vi?2 



So) 



sin(7rA„b) 



d6M{6) 



r(f-A„-n)r(f + A, + n-/i) 



sinf— ttA 
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d6M{6) 



r(f-A. + i)r(f-A, + i) 

r(f - A,-n)r(f + Afe + n-/i)" 

r(|-A, + i)r(f-A, + i) 
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r(f-A. + i)r(f-A, + i) 



f82) 



The two conformal blocks have combined to give the integral in brackets times a real coefficient! 

The case where A^ = A^ can be considered as a limiting case of A^ 7^ A;,, and so is included 
in our proof above. However, one can still ask how the imaginary piece cancels technically 
in this case, since there are no longer two different conformal blocks from the left and right 
side of the diagram to cancel against each other. The resolution of this issue is that this is 



exactly the situation where the derivative relation equation (44) is satisfied. Therefore, each 
conformal block contributes exactly as a total derivative: 

d 



nl ^ ^ 



(83) 



The imaginary piece of the conformal block coefficient eq. ( 73 ) is smooth as an analytic function 
in n, and so when we take the fiat-space limit the sum becomes an integral over a total 
derivative, and therefore it vanishes. 



this characterization of them is completely unaffected by our formal relabeling of the fields. This proves the 
claim. Furthermore, since (f>l s never appear as internal lines, there will be no poles corresponding to them 
in the Mellin amplitude itself - all their poles appear solely in the F functions in the definition of the Mellin 
integrand. This indicates that the appropriate non-perturbative definition of edge cuts is contributions to 
correlators from the poles in these F functions. 
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Next we need to explain why the only contributions to equation (77) are from /c-trace op- 
erators dual to states composed of stable bulk particles. As we mentioned in the introduction, 
although unstable particles can be included in the perturbative cutting rules, the S-Matrix is 
only a well-defined unitary transformation between states made up of exactly stable particles. 
Thus the operator dual to an unstable particle x AdS will not appear in the optical 
theorem, and it must make a vanishing contribution to the conformal block coefficients in the 
flat space limit of AdS. This follows because will mix very quickly with the multi-trace 
operators into which it can decay; roughly speaking, if x has a lifetime r, then it will only 
exist for a time of order t/R, so its contribution in the flat space limit will go to zero. We 
saw this effect in the concrete example of a bulk fi(f)'^x theory in [H], where we derived the 
Breit-Wigner resonance behavior from a re-summed Mellin amplitude. In that case, for any 
flnite AdS scale R the single x mode gave a flnite contribution, but itself became negligible 
as -R — )■ oo, as it was replaced by the continuum of 20 states. Another familiar manifestation 
of this fact is that the delta function resonance from a stable particle is inflnitely sharper, 
and therefore inflnitely taller, than the smooth resonance from an unstable particle. So while 
dropping the single mode corresponding to a stable particle would completely erase its delta 
function resonance, dropping an unstable particle mode has a negligible effect on the S-Matrix. 
In summary: only operators dual to stable particles make an appearance in equation (77). 

Finally, we will complete the argument by showing that the sum over fc-trace operators in 
equations (76) and (77) becomes a. d + 1 dimensional /c-particle phase space integral. 



3.3 Sums Over Operators as Integrals Over Phase Space 

To complete our derivation of the optical theorem and of Cutkosky's 'cutting rules', we need to 
show that the sum over the exchange of all fc-trace operators turns into a phase space integral 
over fc-particle states in the flat space limit of AdS/CFT, as pictured in flgurejsj 

We can understand this by noting that in the large limit, the space of states created by 
single-trace CFT operators is isomorphic to the Fock space of free particle states in AdS. This 
follows from the fact that in AdS/CFT, the Hilbert spaces of the two theories are identical. 
For example, if we quantize a free scalar fleld in AdS [271 Ell EHl [15] , we flnd 

(f){t, p,n) = (f)ni{t, p)anij + 4>*nu{t, fi, p)a\ij (84) 

while the dual CFT operator can be quantized in terms of the same creation and annihilation 
operators, aj^^j and anij-, as 

Oit, n) = Y,^ (e^^'-*y,,(^)an/j + e-^-'*F;,(^])4,,) (85) 

We gave the explicit wavefunctions and normalizations in [15], but the crucial point is physical 
and independent of the details. As is well known, particles in AdS behave as though they are in 
an IR-regulating cavity with a size set by the AdS length R. Thus for flnite i?, the spectrum 
of fc-trace states behaves like the discrete spectrum of fc-particle states in a box of size R. 
When we take the flat space limit R ^ oo, the discrete modes approach a continuum, and we 
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recover the usual fc-particle Lorentz invariant phase space when we sum over these modes. An 
exphcit analysis of the wave functions confirms this intuition [15], and the standard AdS 
quantization above reduces to the flat space quantization of a free field in spherical coordinates. 
Since from AdS/CFT we know that the hilbert spaces of AdS particles and CFT states are 
identical, we can conclude that by summing over a complete set of particle states we are also 
summing over all possible CFT operators. 

As a concrete example, in [16] one of us considered this process in detail for the case of 
double trace operators and 2-particle states. There it was shown that the state created by a 
double-trace operator [O1O2] n/,j can be expressed as 

\n,£,J) = [ dpYuiP) [ dUf{q)\q + p)\q-p) (n » 1) (86) 



{2nYV2RE 

where the state \k) is a one-particle state with momentum k, and the labels J denote various 
angular momentum quantum numbers. The important point is that for primary double-trace 
operators, the function /(g) is fixed to be a Gaussian with width ^JeJR. This means that in 
the fiat space limit where n = ER we obtain precisely the ith partial wave corresponding to 
a 2-particle state with center of mass energy E. 

As we have discussed above, the coefficient of the conformal blocks -Ba,^ at a dimension A 
can be computed by summing over the squares of appropriately normalized 3-pt functions 

{OrO^OM = J2 {cnuYB^,, where cMf,! = (01(^1)02(^2)! n, £, J) (87) 



If we write the double-trace states \n,i,J) using equation (86), then in the fiat space limit 
we can sum over the angular momentum quantum numbers i, J while fixing the dimension 
n ^ ER at the center of mass energy of the scattering process, as measured in AdS units. This 
returns us from the spherical to the plane wave basis, giving 



2 



^ [ ^2nr^ ) J dpdp'Y.YumuipyiniP,-p)cZip\-p) 
^ d k„ d kh r(i+i 



(2T:Y2\ka\ (27r)^2|fcfe| 



where the new 3-pt function coefficient c^^j^{ka, kh) in the plane wave basis can be re-interpreted 
as the square of the flat space scattering amplitude for 12 — ab, and PcoM is the d + 1 
dimensional center of mass momentum. To complete this re-interpretation, the external states 
must also be plane wave scattering states. Plane waves with energy u and velocity v are 
created by acting with [HI IH] 



, d-3 



I"-*) = (2.)^..LW)^- X/«^-'O(*.-f.)|0> (89, 

Thus c^Euika, kh) can only be interpreted as a scattering amplitude when the operators Oi and 
O2 are normalized and integrated in this way. 
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Figure 6: This figure shows how the one-loop example studied in section 3.4 can be computed 
by using the Kallen-Lehmann representation in the bulk of AdS, as discussed in [S]. In this 
way one can write the product of propagators in the loop as a sum over tree-level exchanges 
with different bulk masses, or CFT dimensions. The conformal block decomposition of this 
Witten diagram was indicated in figure |4j 



We should emphasize that the emergence of a ci + 1 dimensional phase space integral from 
the sum over operators is essentially kinematic. It follows as a consequence of the structure of 
the conformal algebra as it reduces to the Poincare algebra in the fiat space limit of AdS /CFT, 
as we discussed in |15]. In particular, this means that if we instead study superconformal field 
theories, then we will instead be taking the fiat space limit of the superconformal algebra. Since 
the spacetime geometry follows from the algebra, in the future one should be able to obtain 
higher dimensional phase space integrals corresponding to decompactifying bulk dimensions 
from the fiat space limit of superconformal theories. 



3.4 A Complete One-Loop Example 

Now let us put the pieces together and understand how the optical theorem applies to the 
one-loop amplitude for 20 — >■ 2-?/' in a theory with couplings f ^iid f x^V^^; as pictured in 
figure |6} We could just as easily treat (j)'^ theory, but we have introduced three fields in order 
to separate out the various different contributions to the one-loop amplitude. 

We computed the relevant one-loop Mellin amplitude in [5] and verified that it has the 
correct fiat space limit. We also explained how branch cuts arise from the coalescence of poles 
in the fiat space limit, and verified that the discontinuity across the branch cut is correctly 
reproduced. Let us now summarize the method and results. To compute a certain class of 
loop amplitudes, one can use the fact that in position space in the bulk of AdS 

GA,(X,r)GA,(X,F) = 5^aA„A2HG'Ai+A2+2n(X,F), where (90) 

n 

^ {h)n (Ai + A2 + 2n)i,fe(Ai + A2 + n - 2/i + 1)„ 
aAiA^l^J 27r%! (Ai + n)i_,(A2 + n)i_,(Ai + A2 + n - ' 

This allows a Kallen-Lehmann type representation [30j for the loop amplitude as a sum over tree 
level exchanges with dimensions 2A^ + 2n, as indicated in figurelo] Taking N/\^{n) = aA^^A^{n) 
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and using the diagrammatic rules from [Qj we find the Melhn amphtude for this diagram is 



M'-'°°'{S„) = At, y;]ViJn)Afci^+2„(iij) (91) 



n 



where M/^{Sij) is the Melhn amplitude for a tree level exchange of a bulk field dual to an 
operator of dimension A. It is 

where 5 = 2A(^ — 2b\2 and A = 2A^ + 2n in our case. The residue Rm is [6l [TJ [9] 

o ^ 1 r(A,-/^ + f)r(A^-/. + f) (i-A^ + f)^(i-Av. + f)^ 

(47r'')3r(A<^-/i + l)2r(A^-/i + l)2 m!r(A - /i + 1 + m) ^ ^ 

and it can be easily computed using the diagrammatic rules from [3 El El [ID] . 

The particular decomposition in equation (91) is very useful for several reasons. First of 
all, to compute the fiat space limit of the loop amplitude we need only know the fiat space 
limit of MA(5jj). But this is simply the fiat space scattering amplitude corresponding to the 
tree-level exchange of a particle of mass (2A^ + 2n)/R. In this limit, the sum over n becomes 
an integral and we can take the large n limit of (n) to find the loop amplitude 

A^i-ioop/ N ^ \g f (^) where Na in) ^ ^ n^^^'^^ f94) 

Another useful feature of equation (91) is that it can be immediately related to the conformal 
block decomposition. As we discussed in |9], conformal blocks and tree-level AdS exchanges 
have identical poles and residues in the Mellin amplitude; they differ only in their asymptotic 
behavior at large 5ij. This means that we can immediately read off the coefficients in the 
conformal block decomposition of the correlator corresponding to this 1-loop Mellin amplitude 
in the [C'x^xl^.o channel, it is 

+2n,o(%) (95) 

n 

where 

P2A,+2n = \gN^^{n)U^~^H-''--''-A (96) 



/(f)'r(A^-f)r(A,-f)r(A, + f-/.)r(A, + f-/.) 
r (A^ + 1 - hf r (A^ + 1 - hf r (^) r (A + 1 - /i) 

with A = 2A-^-|-2n, and we have taken £ = because the coefficient of all the blocks with £ > 
vanish. The normalization here follows from the relative definition of the Mellin amplitude for 
a bulk exchange and the normalization of the conformal blocks, which is most easily determined 
by relating the residues at their poles. Finally, one can also see [9] from equation (94) that the 
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discontinuity across the branch cut in the flat space amphtude is just given by the residue of 
the pole in this equation, which is 



disc [M^- 



N, 



loop] 



Xg- 



(97) 



This formula appears on the left-hand side of the optical theorem, as the imaginary part of 
the scattering amplitude. It remains to see how this is reproduced by the right hand side of 
the optical theorem. 

To compute the right-hand side of the optical theorem and check equality, we need to 
apply our conglomeration procedure to compute the OPE coefficients that determine iVA^(n). 



2.3.4 



and found the OPE coefficient 



Fortunately, we did this calculation in section 
for the operator [O^(9^]„_o in the product x C^. One can verify using equation (60) that 



6c. 



6ct 



P, 



(9J 



using the expression we computed in equation (96). This provides a very non-trivial test of 



the conglomeration technique. We can also interpret the product of OPE coefficients in terms 
of a 2-particle phase space integral in the fiat space limit. By definition, we have that 



(/)(/) rpn,£,J 



(0^(Pi)0^(P2)|2A^ + 2n,£, J) 



(99) 



If we apply equation (89) to the operators, then we see that at large n, the constant c^^j 
can be interpreted as a flat space scattering amplitude between two particles in plane wave 
states and a 2% particle state in the spherical wave state of equation (86). So each OPE 



coefficient can be interpreted as a tree-level scattering amplitude for 20 — )■ 2% and 2il) — 2%, 



respectively, and in the flat space limit equation (98) becomes the optical theorem. 



4 Discussion 

Although it has been clear for some time that the S-Matrix is the only exact observable in flat 
space quantum gravity, it has remained somewhat mysterious what sort of holographic theory 
pn 123] might compute it. A true theory should do more than just output scattering 
amplitudes, it must also provide an understanding of how fundamental principles such as 
bulk locality and quantum mechanics emerge. This would appear to be an especially difficult 
problem if one tries to obtain a holographic theory that 'lives' directly on the null boundaries 
of flat spacetime, where notions such as time and distance can only have a limited meaning. 

We have argued that the flat spacetime limit of AdS/CFT may provide the theory that we 
have been looking for. In [9] we derived a formula conjectured by Penedones [6] that relates 
the Mellin amplitude [H [5] for n-pt CFT correlators to the flat space S-Matrix. This formula 
expresses the S-Matrix as a simple integral transform of the Mellin amplitude, which itself must 
be a meromorphic function restricted to have only simple poles on the real axis. The clearest 
way to understand flat spacetime locality from a holographic perspective is via the analyticity 
properties of the S-Matrix. This strongly suggests that the most natural way to understand 
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locality may be in terms of the very restricted analytic structure of the Mellin amplitude 
combined with some set of assumptions about the spectrum of the CFT, as described in [22j . 

In the present work we have shown that the operator product expansion and conformal 
block decomposition of CFTs encodes unitarity in a form that appears very similar to the 
usual optical theorem for the S-Matrix, as pictured in figure |2} By taking the flat space 
limit of the AdS/CFT duality we showed that one can derive the optical theorem directly 
from these unitarity relations for the case of 2-to-2 scattering of massless scalar particles. In 
the perturbative case this reduces to the usual cutting rules for Feynman diagrams. We also 
saw something subtle and interesting occur with the spectrum - in the fiat space limit, CFT 
operators dual to unstable particles must decouple from the unitarity relation. It would be 
interesting to understand this phenomenon better, and to explain it without having to appeal 
to bulk reasoning. To give a vaguer and more ambitious-sounding summary, one might say 
that we have shown how to derive bulk quantum mechanics from the quantum mechanics of 
the holographic dual. 

To make these derivations possible we developed technology to conglomerate many local 
operators together into a single composite operator. To enact this conglomeration we used 
smearing functions /a/ that we labeled 'wavef unctions' because, via the operator-state corre- 
sopndence, these functions also extract definite states from the CFT. The relative simplicity 
of our formalism was made possible through use of the Mellin amplitude for CFT correlators. 
Since the Mellin representation depends on position space kinematics through specific power- 
laws, it accords naturally with the structure of the conglomerating wavef unctions, which have 
the power-law behavior of CFT 3-pt correlators. Perhaps in the future the logic will be re- 
versed; one might attempt to derive the Mellin amplitude as the representation where these 
wavefunctions behave most naturally. 

We also gave a schematic argument that the usual Dyson series for the S-Matrix can be 
constructed in AdS using the dilatation operator, and that at first order it gives rise to the 
prescription for the S-Matrix in terms of the anomalous dimension matrix that was given in |16j . 
In [9] we derived Penedones formula for the S-Matrix using essentially the same wavepacket 
setup that was originally described in [121 US] and further examined in [HI [15] . Thus we have 
given at least a rough explanation why all three of these holographic formalisms for computing 
the bulk S-Matrix agree, although we view the formula in terms of the Mellin amplitude as by 
far the most elegant and physical. 

The impediment to formulating a general proof of the optical theorem for any n-pt scatter- 
ing amplitude with particles of arbitrary spin seems to be mostly technical. The requisite CFT 
technology to describe higher-spin particles [25l [20] and higher-point conformal blocks has not 
been fully developed, so we lack the necessary CFT ingredients. Note that the fact that sums 
over /c-trace operators reduce to fc-particle phase space integrals is essentially kinematic, fol- 
lowing from the structure of the conformal algebra when applied to large dimension operators. 
Thus in the future one should be able to show how superconformal theories give rise to higher 
dimensional phase space integrals corresponding to extra dimensions that decompactify in the 
fiat space limit. 

Perhaps it would be appropriate to mention here that some of the standard lore concerning 
n-pt correlators seems a bit misleading, and that this becomes very apparent when it is re- 
stated in terms of the bulk S-Matrix. CFTs are often viewed as extremely constrained theories. 
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and the statement is often made that all CFT correlators are completely determined by the 

3- pt correlators. In 2-dimensions where there is an infinite dimensional symmetry algebra this 
may be a very powerful point, but in higher dimensions it is rather trivial. In particular, the 
flat space limit of this statement is the trivial claim that once we know the 2-to-fc scattering 
amplitudes for all k, we know the entire S-Matrix. 

We have mentioned this point in order to emphasize that higher correlation functions in 
CFTs are very non-trivial and are worthy of investigation. Higher-point conformal blocks have 
rarely been discussed, but we would require such objects in order to formulate S-Matrix uni- 
tarity beyond four particles. It would be interesting to develop techniques for handling them. 
Another reason to study higher-point correlators is to understand Hawking radiation, which 
produces large multiplicity final states that can be understood in terms of the properties of 
correlators involving a large number of single-trace operators. The fiat space limit of AdS / CFT 
turns the bootstrap program for CFTs into the S-Matrix program, so focusing solely on the 

4- pt correlators of single-trace operators would be like restricting the study of scattering am- 
plitudes to 2-to-2 processes. We hope that our portrait of the holographic S-Matrix suggests 
new directions for future investigation. 

The AdS/CFT correspondence [H Ej [3] has proved to be an extremely important and 
fruitful discovery. One perspective on this correspondence views it as a very large compendium 
of exact dualities between various CFTs and other models of quantum gravity, often in the 
form of Superstring Theories or M-Theory. A complementary perspective involves placing an 
arbitrary effective field theory in Anti-deSitter space and computing the correlation functions of 
bulk fields as they approach the boundary [161 |3l]. The correlators of this 'effective conformal 
theory' wfll approximate those of a CFT to very good accuracy [22], with errors suppressed 
by powers of operator dimensions divided by the cutoff in AdS units [16]. The search for a 
non-perturbative completion to a gravitational EFT in AdS space translates into the question 
of whether there exists an exactly defined CFT that approximately reproduces the boundary 
correlators of the AdS effective theory. Similarly, we can UV complete a gravitational EFT in 
fiat spacetime if we can find a sequence of CFTs of increasing central charge whose spectrum 
and correlators approximate those of the EFT in the fiat space limit of AdS. 

In other words, although it may be extremely challenging to actually find a non-perturbative 
completion for a given gravitational effective field theory, anyone can use AdS/CFT to correctly 
formulate the question. This means that one can obtain robust results about quantum gravity 
by modeling the correlators of low-dimension operators using bulk effective field theory, and 
then using the bootstrap approach [TTl [191 ESI |36l |25l [201 [33 [Ml [391 126] to constrain correlators 
that involve operators of larger dimension. Since dimensions in the CFT correspond to bulk 
energies, one can obtain information about processes at trans-Planckian energies in AdS. In 
[9] we used Hawking evaporation to make a prediction for the conformal block decomposition 
of 4-pt correlators. If one can derive this generic behavior from CFT dynamics by using the 
bootstrap, then quantum gravity may be accessible to mortals. 

Should we view a holographic description as the final word on quantum gravity in a par- 
ticular class of spacetimes? The legalistic answer may be yes, but it seems that holographic 
descriptions such as AdS/CFT do not readily yield information about the physics behind hori- 
zons, and we might hope that such questions are not entirely ill-defined. It seems reasonable 
to assume that if we could do experiments on black holes in a large enough laboratory, we 
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would see unitary evaporation, with intrinsic errors that decrease as we increase the size of 
our detector. This suggests that it may be worth looking [IQ] for an approximation scheme 
beyond effective field theory that encodes both locality and its demise. 
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A Conglomerating Operators: Regular izat ion Details 

Here, we will describe in more detail our regularization procedure for conglomerating operators. 



Let us return to eq. (39), but keeping A arbitrary: 



_ .. r(^^)r {h - AO r {h - a,) r {-h + ^±^) 
ca.o^Vo ^ r(Ai)r(A2)r(/i + ^=^^)r(2/i-^±^) ^ ^ 

To fix the normalization, we want to conglomerate operators again to obtain a two-point 
function. As we will see in a moment, it is necessary to regulate by both taking the dimension 
A' of the second operator in the two-point function to be arbitrary and also to define it as a 
function of two positions P5 and P5/ that we will take to be equal in physical quantities. That 
is, we will calculate the two-point function 

{OA,o{Po)OA'fi{P5,P5')) (101) 

where 

0A',o(P5,P5') = T^ [ d'P^d'P2 .,.. -r^ Oi{Pi)02{P2) (102) 

^V.O-^ A2 ' Pi, ' ^25' ' 



Then, combining (39) and (102), a straightforward application of Symanzik's star formula 
demonstrates that 



(0A,o(Po)0A',o(n,P5')) = Jd''P,d''P2 



A' + Ai+A2-2d 

P12 ^ I CaiCa2-Pi2 



Ai+A'-A2 A2+A'-Ai I A+A1-A2 A+A2-A1 
02 02 \ p 2 p 2 

^15 ^25' \-^01 ^02 



CaiCa2 



p 2 p 2 p 2 
-'05' -"OS 55' 

Ca,0 Vr2'^r(^)r(/l - A+A2-Ai )p(^ _ A+^^)r(^ _ 

""ivj^ r(/i)r( ^+^]-^^ -)r( ^+^r^0 r(2/i- A) ' 

(103) 
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where we have taken the hmit A' — )■ A in places where it does not produce singularities. Now, 
let us calculate the inner product of the states corresponding to Ca,o(-Po) and Ca',o(-P5, -Ps') 
in the usual way in radial quantization by taking xq — )■ and x^^x'^ — )■ oo and rescaling by 
the appropriate powers. As we take A' — > A, this inner product has the interpretation of a 



normalization, so this fixes A^^q (^l)? but with "0" replaced by ^' ^ ^ : 



r(/i)r(^±^)r(^±^)r(ci-A) ' ^ ^ 



A' = Ai + A2 + 2n and A = ^7 + A2 



Finally, putting together eq. (100) and (104) to obtain the physical quantity c^qj ^ay take 



2n + e and safely take the limit e — )■ 0, for which we 



obtain the finite result of eq. (42) 



B Conglomerating Operators: Spinning Conformal Blocks 

We will extend the calculation in section [2] to general spin. The general form of three- and 
two-point functions of primaries is 



-'12 

mP,)0,{P,)[OsO,U+e{Ps)) = gf^^ ^^"';,1'!.a1~Xa" f£f • (105) 



,{{Z, -Pi)P2S-{Z^- P2)Pl3 

3-n pA2+n+^pA 
12 -'23 -'31 

The general spin smearing functions are again just three-point functions with shadow fields: 



[O:OJ„,,(/>0) = A /<i/'i<i/^.p^;^5^^i^^^fS0i(A)0.(/>.). (106) 

■J -'12 -'02 -'01 



We will be interested in conglomerating Oi and O2 in order to obtain spinning operators in 
their OPE. While this procedure works for extracting a general operator, it will be interesting 
to apply this to special CclSGS cLS well. In particular, our first application will be to the theory 
at infinite N, where the only operators that arise in the OPE are the double-trace operators. 

B.l Conglomerating Operators: Disconnected Four-point Function 

Since the only operators to consider in this case are the double-trace operators, we will label 
them by their indices n, i, where A = Ai -|- A2 -|- 2n -f i.e. we will use the notation [Oi02]n,e 
rather than [Oi02]a/- It will also be convenient to define A^ = ^J^^. To extract the N = 00 
three-point function of [Oi02]n/ with O2, we integrate the smearing function against the 
disconnected diagram: 



{[O,O2]nAP0)OliP3)O2iP,)) = ^ J dPldP2 



{{Zo ■ Pi)P2o - (Zo ■ P2)PioY CaA 



pd-2Aa-n pAi+n+e pA2+n+e pAi pA2 
n,£ ^ -'12 02 -"01 13 -'24 



(107) 
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To do this integral directly is more complicated than the ones we have encountered, because 
the contraction vector for the indices of the spinning field acts like an additional point. To 
perform this integral more simply we can use the fact [251 120] that the three-point function of 
scalars with a spin-£ field has to be made out of powers of the tensor 

Cqab = ZqPq - Z^P^, (108) 

and the correlator can be written J^tp°oAp'--iA ■ particular, the correlation functions have the 
gauge symmetry F(Pj, Zi + aiPi) = F{Pi, Zi). Thus, one way to make this constraint manifest 
and simplify the calculation is just to pick a gaugej^ Of course, we have to decide what gauge 
to choose. Seemingly natural choices like Z^ ■ P^ and Z^ ■ Z^ are not helpful, since they are 
already satisfied in any gauge. The next best thing seems to be Z^ ■ Pj for one of the other j's. 

How much does this simplify the computation of the OPE coefficient? Let us take the 
gauge where Zq ■ P4 vanishes. Now, any time we get contributions from a positive power of 
Zq ■ P4, we can drop them. So far, this does not use any knowledge of the final form of the 
correlation function other than the fact that it is gauge-invariant; for any Zq, we can always 
re-obtain the full gauge-invariant result by calculating in this gauge and then restoring gauge 
invariance by taking ^0—7-^0 + oPq, a = which is simply the gauge transformation 

that took Zq ■ Pi = 0, and which clearly reintroduces the gauge-invariant P3 ■ Cq ■ P4. 

^ In order to gauge- fix the Zi^s, we should first decide what these abstract objects look like. They are clearly 
not regular points like P^'s in the boundary theory, because they satisfy Zi - Pi = 0, but Zi ^ Pi (otherwise 
the C/^^'s would vanish). Let us imagine for a moment however what they would look like if they were, so Zi 
projects down to the point z^. Then, Zi ■ Pi = implies 

iz,-x,f = 0. (109) 

Obviously, in Euclidean space this requires Zi = Xi, which we do not want. Without loss of generality, let's 
take Xi = 0. Then, we need 

zf = 0. (110) 

Now we see that Zi should simply be a point on the boundary with complexified coordinates. So, in general, 
we need 

Zi = Xi + qi, qi=0- (111) 

This also makes it very explicit why taking Zi to be one of the Pj 's in order to simplify the integrals is not 
allowed. We can work out how the gauge transformation Z — >■ Z + aP acts on q: 

Z = (1, z^, z^) ^{l,x-ix + 2q),x^' + q^) ^^-^"^ {l + a,x- ((1 + a)x + 2q), (1 + a)x'' + q^') 

= {l,x-{x + 2-^),x'^ + -^). (112) 
1 + a 1 + a 

In the last line, we have used the fact that P's project down to boundary points by rescaling P+ — >■ 1. Thus, 
the effect of the gauge transformation on Z is written in terms of q very simply: 

Z^Z + aP ^ q^q' = —^. (113) 

1 + a 
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Having fixed a gauge, we can continue with the computation. The smearing integral to 
perform is 



where Pi^ — — 2Pj • Zq. Focusing on individual terms, we compute 



/pfc pi~k 



Bek — / ^-Pl'^-P2 ^_2A„-npA2+n+fcpAi+n+f-fcpAipA2 
-'12 -'01 -'02 13 -'24 



( Ai + n + ^ - k)s,^ {d - 2Aa - n)s,^ ( Aa)^,, 

/(Ai + n + £ - A; + 5o^, d - 2A„ - n + 5^, A2 + 54^), (115) 



where we have introduced the notation 



I{a, b,c) = [ d'P, \ = 6, c)i?,"4,o<,i^oi,4, = f ' > 

^^'''''''^^ - r(a)r(6)r(c) " ^ 

The ^ix's satisfy the constraint 5qz + (^iz + (^4z = — A;, and the arguments of the / function get 
shifted since we have eliminated 3 5y variables by using the constraints, 5o4 + ^14 = A2 + ^2^, 
etc. Now, the only poles that survive our gauge choice are 5qz — ^iz — 0, so — — k, and 
we have 



Bek — 



dPiPlxi^i +n + i-k,d-2A^-n + £-k,A2){d-2Aa- n)e-k 



pAi ph+n+e pAi+A2+n-h ph-Ai-n 
-'13 -'01 -'04 -'14 



x(Ai + n + £- k,d-2Aa-n + i- k, A2)x{h + n + £, /i - Ai - n, Ai + £) ^ 



pA2+n pAi+n+^ p-n 3z 
-^04 -^03 -^34 

x{d - 2Aa - n)e-k{Ai)i 

(117) 

Somewhat remarkably, performing the sum over k obtains a relatively simple result for the 
three-point function: 

mP,)0,{P,)[0,0,UUPo)) = cll 



pA2+n pAi+n+l* p—n ' 
-^04 -"03 -"34 



_i2 ^ CaiCa2 7r^^r(-n)r(/t - Ai)r(/i - A2)r(-/i + 2A„ + n + £) 
~ 7V„^^, r(/i + n + £)r(Ai)r(A2)r(2/i-2A„-n) 

(118) 



This indeed is the correct form of the three-point function in • P4 = gauge. 
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Having obtained the three-point function, we next need to conglomerate again in order 
to determine the two-point function and thus the normahzation of [Oi02\2n+e- The smearing 
integral we have to compute is 



'^ril f IT-, IT-, {^5 ' P3P45' ~ ■ P4P35) 



r,d-2Aa-n' pA2+n'+£ pAi+n'+^ 
n/ ^34 ^35 ^45' 

{Zq ■ P3P04 - Zq ■ P4P03Y /,,pv 

^ pA2+n+^pAi+n+£p-n ' \^^^) 
-^04 -^03 -^34 

where we have done the usual n ^ n',5 5' regularizations. We will take the gauge Po • Z5 — 
0,Pr, ■ Zq = P5/ ■ Zq = (the second two are the same because P5 = P5/ everywhere except 
when they may give rise to P55'). Binomially expanding, we have 



{[O,O2]2n+e{P0)[O,O2]2n'+e{P,)) = f"^) E ( ^ ( 



Be.kk' = / dPzdPi 



n/ k k 

TDk' nt-k' pk pi-k 

3^'-* 4^' S^-* 4z 

pd-2Aa-n'-n pA2+ri'+fc' pAi+n'+£-fe' pA2+n+^-fe pAi+n+fc ' 
-'34 35 -'45' -^04 -^^03 

(120) 

We first perform the P4 integration, which results in the introduction of five new d integration 
variables after imposing constraints. Four of them are forced to vanish by our gauge choice, 
leaving a single variable that can be converted by the residue theorem into a sum: 



ekk' 



^ V 1 f .p Pi;rp!,:"^PT, , ..^+1 {i-kw-h')\ 

p-h+2A„+n+„'+m j 3 A2+n'+fc' pAi+n+fc I ) rn\{i - k - my.{i - k' 



X 



05, J ^35 ^03 m\{£ - k - m)\{£ - k' - m)\ 

X{d - 2Aa -n' -n + 2i-k-k' -2m,Ai + n' + e-k',A2 + n + e-k) 

ph— Ai—n' -\-i—k—m ph— A.2—n-\-l—k' —m 
-"03 35' 

X (d - 2Aa -n' - n)2i-k-k'-2m- (121) 

We can simphfy this by taking P5/ — )■ P5 in places where it will not lead to singularities, which 

in particular is any place that does not have a n or n' exponent. So, we can exchange the 
powers of P35 and P35/ in the denominator for Pj^^"' ^^^^Pg^". But, then there are no powers 
of k, k' remaining in the Pj^'s, so we can complete the sum over them outside the integral. We 
can furthermore take n — > n' in the prefactor, since this is needed as a regulator only in the 
powers of P^^'s. We thus obtain 



{[OlO2]2n+e{P0)[O,O2] 



" 7i''{-l)"'+\i\)^T{h-n-A,)T{h-n-A2)T\i + 2n + 2Aa-h) 



-2Y^N\ m!r2(l + m)r(Ai + (i + n)V{A2 + £ + n)T{2h - 2n - 2A„)r(m + 2n + 2A„ - h) 

/pi—m p£—m pm 
fjp ^3z ^3z' ^zz' /-|22^ 
^ ph+n-n' -\-t-m p-n ph-\-n' -\-i-m p-/i+2Aa+n+n'+m ' V ' 

^03 35' 35 ^05' 
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This last line is exactly of the form of Symanzik's star formula: 



r)t-m T}t-m pm 

^ ^3z ^3z' ^zz' 

^ ph+n—n'+£—m p— n ph+n' +l—m p—h+2Ha+n+n' +m 



03 35' -"35 05' 

_lY+m+i pi_^^(^^ _ ^)!^(/^ ^n' + l-m,h + n-n' + ^-m,-n) 

2/^a+2n'+l pn'-n ' 



(123) 



pziSa+'!.n' +1 pn- 
-'05 -'55' 

where we have taken P^i — )■ P5 except inside P55' as well as n — )■ in some places that do not 
lead to singularities. Putting everything together, we thus have 

2£ 



pi 



^\2J p2Aa+2n+£' 



cly TT^^T{0){-iyi\T{h - n - Ai)r(/i -n- A2)T{21 + 2n + 2A„ - l)T{-h + I + 2n + 2AJ 
ivjr Y{h + l)T{l + n + Ai)r(/ + n + A2)T{2h - 2n - 2Ajr(/ + 2n + 2A„ - 1) ' 

(124) 



From equations (118) and (124), we can choose the normalization factor Nip to set a canon- 



ically normalized two-point function coefficient, 03'^ = 1. Then, we can read off the OPE 
coefficient, which after some simplification is 

._12x2 ^ i-lYCA,CA,{Ai -h+ 1)„(A2 -h + l)„(Ai),+„(A2)^+n 

^"""''^ i\n\{i + h)^{A, + A2 + n - 2/i + l),(Ai + A2 + 2n + £ - l)z(Ai + A2 + n + £ - /i)„' 

(125) 

B.2 Conglomerating Operators: Connected Four-point Function 

Next, let us apply the conglomerating methods to connected four-point functions. Here, it is 
more convenient to work with the Mellin representation of the four-point function, since all 
position space integrations can be performed using Symanzik's start formula. We will use as 
our Mellin coordinates the variables 

6 7 + (5 

x = A„-5i2 = -, y = Su = ^—, (126) 

since our expressions will typically be more compact in terms of these variables than in terms 
of 6i/s or S,j. The four-point function can then be written as 

M{x,y) 

r ' ^ ^ 1 

J ^ ^ '^^ ^ " ^ ^ ' ^ ^ ^' Pt2-''Ptr^{Pi3P2,Y-y{PMy 

(127) 

where M{x,y) is the reduced Mellin amplitude. We obtain the OPE coefficient by smearing 
to produce a three-point function: 

= / dPidP2 ^ p °-2Ai pA„in+° pA!ln+] (128) 
J -'12 -'02 -'01 
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We will again work in a gauge where ■ P4 = in order to simplify the calculation. In order 
to use Symanzik's formula, we will binomially expand to obtain 

B = [ dxdvdP dP M(x, y)PizP2z q29) 

I r>Aa+n+k pAa+n+e-k pd-Aa-n-x pAt-x / p p W-y/p p \y ^ ' 

J ^02 ^12 ^34 lA3-f24j ^l-ri4-f23j^ 

Performing the Pi integration first, using Symanzik's formula, we obtain 



m=0 



Bimk = [dS]dxdydP2M{x,y) 



pi—m pm 
^2z -^3z 



pSo2+^a+n+i-k p(5o3 pSo4 p&23+y p&2i+x-y pSzi+^b-x 
-'02 -'03 -'04 -'23 -'24 -'34 



r(^02)r(^03)r(^04)r(^23)r(^24)r(^34) .^3^. 

V{K + n + k)V{d- K-n-x)V{x-y)V{yy ^ ' 

The P2 integration also follows from the application of Symanzik's formula. The essential 
structure of the result of conglomerating is therefore that it introduces a projection function 
i7(A,£;x,y): 

^3 = pA.+n+£pA^+npA,-A.^n / dxdyM{x,y)V{/\,-x)V{/\,-x)H{A^ + n,l-x,y), 

-'03 -'04 -'34 J 

(131) 

where an important point is that H{A, i; x,y) does not depend on or A;,: 

f . . r(^04)r(^34)r(A - 6o^)Tiy - <5o4 - ^34)r(x - A - ^34)r(/i + k-m-x + ^34) 

°^ '^^T{k + A)T{x-So4-S34)T{d-A-x)T{h-A + 6o4)T{h + l-m + A-x + 634) 
xT{h - A-y + 6o4)T{h + l-m-x + 6o4 + 634)T{-h + m + A + x-6o4- 634). (132) 

We do not have a nice closed-form expression for H{A,i;x,y) in general like we do for the 
special case of / dyr{y)r{x — y)H{A,0] x,y) that appeared in section [2j It is possible that 
such an expression exists and could be obtained with more effort, and could be useful for 
extracting OPE coefficients for specific theories. In addition, one may perform the 6q4 and ^34 
integrations above using the residue theorem in order to obtain H{A,i;x,y) as a sum, which 
could perhaps be useful in some situations for numeric computations. 
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C Double Trace Operators 



The purpose of this appendix is to obtain a completely general recursion relation that expresses 
double trace primary operators in terms of linear combinations of derivatives acting on Oi02- 
The full conformal algebra is 

[M^i,, Pp] = i{r]ppPu - VupPp.), [Mpy, Kp] = i{r]ppKy - r]^pKp), 
[M^,, D] = 0, [Pp, K,] = -2{7]p,D + iM^,), 
[D,Pp]=Pp, [D,Kp] = -Kp. (133) 

Primary operators are those annihilated by the special conformal generator Kp, so that 
[Kp, O] = 0. We will work with operators that are eigenstates of D and the angular mo- 
mentum generators, so our double trace operators can be written as [Oi02]n,e- In what follows 
we will generalize the computations of [21], [6j. 



C.l Action of on T{k, I — k, ui, U2, m) 

By contracting with traceless symmetric polarizations V , we can write the action of Kp on 
the double-trace operators in a particular basis. Let us follow the notation of [6] and write a 
general operator of the desired form as 

T{k, I - k, u„ U2, m) = 'P„, . . . P^^Pp, ...Pp^ (P'rO.P^,^, . . . P^,Pp, . . . Pp„^ {PY'02. 

(134) 

Then, the action of Kp on this operator is as follows: 

KpT{k, I - k, u,, U2, m) = V- [2u^{d - 2u^ - 2A,)PpP^^ . . . P^^Pp^ . . . Pp^{P')^^-'0,{. . . 

( 

-2(Ai -hm-h A;-h2Mi - 1 




+ ^ . . . P^, P^^ ...P^^ ... P^'" 



Ms 



(P^)"iCi(...)C2 



\ s>r s,r s>r / 



+(1 ^2,k^l-k)]. 



(135) 



Here, we have defined P^*\ P^:'^'^ as indicated, and (. . . indicates the O2 half of the double- 
trace operator before the Kp action. Performing the contractions with V and symmetrizing, 
we can write the results in terms of 

Tp{k, I - k, u,, U2, m) = Vp-^^P^, . . . P.,Pp, . . . Pp^ (P'rO.P^,^, . . . P^.Pp, . . . PpjP'rO 



(136) 
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where V"^ is also symmetric and traceless. Grouping like terms, we find 

K^T{k, I — k,ui,U2,m) = 2ui{d — 2ui — 2 Ai) x T^{k + 1,1 — k,ui — l,U2,m) 

2u2{d ~ 2u2 — 2A2) X T^(A;, / — A; + 1, Ml, ^2 — 1, m) 
-2fc(Ai + m + A; + 2ni - 1) x T^{k - 1,1 - k,ui,U2,m) 
-2{l - k){A2 + m + I - k + 2u2 - 1) x T^{k,l - k - l,ui,U2,m) 

— 2m(Ai + m + 2^1 — 1) x T^{k,l — k + l,Ui,U2,m — 1) 
— 2m(A2 + m + 2^2 — 1) x Tfj_{k + 1,1 — k,Ui,U2,m — 1) 

+171(171 — 1) X T{k + 1,1 — k,ui,U2 + 1,171 — 2) 
+m(m — 1) X T{k,l — k + l,Ui + l,U2,i7i — 2). 

(137) 

Since we are demanding that the operator [Oi02]n,e = ^ a{ki, k2, ui, U2, m)T{ki, k2, ui, U2, m), 
be primary, we obtain an equation for the coefficients: 

= 2{ui + l){d-2ui-2-2Ai)a{ki-l,k2,ui + l,U2,ra) 
+2{u2 + l){d-2u2-2- 2A2)a{ki, fca - 1, mi, M2 + 1, m) 
-2{ki + l)(Ai + m + ki + 2ui)a{ki + 1, k2,Ui,U2,m) 
-2(/c2 + 1)(A2 + 171 + k2 + 2u2)a{ki, ^2 + l,Mi,^t2, fn) 
— 2(m + l)(Ai + m + 2ui)a{ki, k2 — 1, Ui,U2, m + 1) 
-2(m + 1)(A2 + m + 2u2)a{ki - 1, k2, Mi, ^2, m + l) 
+ {m + l){m + 2)a{ki - 1, ^2, ui,U2 - l,irL + 2) 

+ {m + l)(m + 2)a(A;i, k2 - 1, Mi - 1, U2, m + 2) (138) 

Note that there are two terms here where the total spin has been incremented to £ + 1, while in 
the remaining terms it has been decremented to £ — 1. These two types of terms must cancel 
amongst themselves. The two incremented terms imply the equation 

= -2(A;i + l)(Ai + m + fci + 2Mi)a(A;i + 1, A;2,Mi,'W2,'"^) 

-2{k2 + 1)(A2 + m + ^2 + 2u2)a{ki, k2 + 1, ui,U2, m). (139) 

This is very constraining, since it completely fixes the fc-dependence of a. Making an ansatz 

a{k,l - k,ui,U2,n - ui - U2) = Sn,i{k)h{ui,U2) (140) 

we can solve for Sn,i{k) uniquely: 

(-1)'= 

^""''^^^ " k\{l - k)\r{Ai +n + ui-U2 + k)r{A2 + n + U2-Ui + l-k)' ^^^^'^ 

Note that this agrees with the results of [U] for the n = case he computed. Substituting this 
back into our constraint on a for the spin i — 1 terms, we obtain an equation for b{ui,U2). 
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After some simplification, this can be written 

Q _ (n - til - U2){1 + n - Ml - U2)h{ui - 1,^2) _^{n -ui- U2){1 + n-ui — U2)h{ui,U2 - 1) 



k{I. - A; + n - 1 - Ml + M2 + A2) {k - l){k + n - I + Ui - U2 + ^ 



2{n - Ml - U2){n - 1 + Ml - M2 + Ai)6(mi,M2) _^ 2(n - Mi - M2)(n - 1 - mi + M2 + A2)&(mi,M2) 



A;(A; + n - 1 + Ml - M2 + Ai) (/ - - + n - 1 - ui + M2 + A2) 

2(m2 + l)i-d + 2(1 + M2 + A2))6(mi,M2 + 1) ^ 2(mi + l){d- 2(1 + mi + Ai))6(mi + 1, M2) 



k{l — k + n — 1 — Ui + U2 + A2) {k — l){k + n — 1 + Ui — U2 + Ai) 



Because of the /c-dependence of s{k), this has many terms that depend exphcitly on k. But, 
we have just proven that 6(mi,M2) cannot have any A;-dependence! Thus, we can muhiply 
through by all the terms in the denominators and collect coefficients by powers of k. This 
gives us exactly three equations for b{ui,U2), one for each of 1, k, k"^ . It is straightforward to 
check that one linear combination of these three equations vanishes, so in fact we obtain only 
two equations. We can take linear combinations of the remaining two equations to obtain two 
recursion relations, one that increments mi and another that increments M2: 

1 N _ m((l + m)6(Mi,M2-l)-2(/ + m-l + 2Mi + Ai)6(Mi,M2)) 
^^''^^^^''^^ - 2(Mi + l)(-rf + 2(l+Mi + Ai)) 

m((l + m)6(Mi-l,M2)-2(/ + m-l + 2M2 + A2)6(Mi,M2)) 

'^^-^^ + '^ - 2(M2 + l)(-rf + 2(l+M2 + A2)) ^'''^ 

where we define m = n — Mi — M2 for concision. This allows us to obtain the full solution for 
a{k, I — k, Ml, M2, m) for any n, I, up to a single overall normalization factor by beginning with 
6(0, 0) and then recursively increasing the parameters mi and M2. 

C.2 Solving in the Boundary Case 

Unlike in the case of Sn,e{k), we are unaware of a full closed form solution for 6(mi, M2). However, 
the equations for b{ui, U2) simplify when M2 = or mi = 0, giving 

^ ^ 2(mi + 1)(-/i+(1 + mi) + Ai) ^ ^ 

and similarly with 1^2. This has the simple solution 

^^^I'UJ ( ^) 2-r(Mi)r(n + l-Mi)(Ai-/^)„, 
when normalized so that 6(0,0) — 1. 

C.3 Large u,n behavior of the coefficients 

The solutions to the above recursion relations are seem to be very complicated in general. 
However, we will look for simplifications at large n. Our first step will be to obtain an recursion 



(142) 
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relation for just the diagonal elements, ui = U2 = u. This may be done by moving along the 
diagonal and near-diagonal ui = M2 — 1, solving only for these elements and no others. Formally, 
we can take h{u, u) = bd{u),b{u — l, u) = bo{u); then, eliminating bo{u) will give us our recursion 
relation for the diagonal elements bd{u). This may be easily done by shifting the argument 
where necessary, and we find 

= Wbd{u){~{2{Ai + u)-d)){2{A2 + u)-d)--^bd{u-2){n-2u + l)i 

+ - 2u) + ^(2Ai + 2l + 2n + 2u- 4)(2A2 + 21 + 2n + 2u - A) - 2{2u - !)(/ + n - 1) + 2u'^^ 

xhd{u - l){n - 2u + 1)2. (146) 

Now, we want to take the large u,n limit of this. At leading order in this limit, the shifts 
in u make no difference and we simply find 

^ n^{3n-8u)bd{u), (147) 

which implies that bd{u) is peaked around u = To go to higher orders, we can expand the 
shifts in u as derivatives. We further take u = ^ + 6, with 6 of 0{1). The subleading behavior 
of the recursion relation is 

^ n3(4(23 + 3c/ + 2/-2(Ai + A2) -32(5)6rf((5) -5n6'^((5)). (148) 

This is a first-order ordinary differential equation, and is easily solved. It has a simple approx- 
imate solution in terms of a Gaussian: 

bd{d) oc exp I — — 

2 5n ^ 23 + 3c/ + 2/-2(Ai + A2) 
^ = 128' = 32 • 

Thus we see that these coefficients are sharply peaked around ui = U2 = ^ for large n. 
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